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Abstract: The holographic dual of a finite-temperature gauge theory with a small number 
of flavours typically contains D-brane probes in a black hole background. We have recently 
shown that these systems undergo a first order phase transition characterised by a 'melting' 
of the mesons. Here we extend our analysis of the thermodynamics of these systems by 
computing their free energy, entropy and energy densities, as well as the speed of sound. We 
also compute the meson spectrum for brane embeddings outside the horizon and find that 
tachyonic modes appear where this phase is expected to be unstable from thermodynamic 
considerations. 
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1. Introduction 

In a broad class of large- N c , strongly coupled gauge theories with a holographic dual, a small 
number of flavours of fundamental matter, N f -C N c , may be described by N { probe Dq-branes 
in the gravitational background of N c black Dp-branes [1]. At a sufficiently high temperature 
T, the background geometry contains a black hole [2]. It was recently shown that these 
systems generally undergo a universal first order phase transition characterised by a change 
in the behaviour of the fundamental matter [3]. 1 

Minkowski embedding Critical embedding Black hole embedding 



Figure 1: Various Dq-brane configurations in a black Dp-brane background with increasing temper- 
ature from left to right. For low temperatures, the probe branes close off smoothly above the horizon. 
For high temperatures, the branes fall through the event horizon. In between, a critical solution exists 
in which the branes just 'touch' the horizon at a point. 

From the viewpoint of the holographic description, the basic physics behind this transition 
is easily understood. Increasing the temperature increases both the radial position and the 
energy density of the event horizon in the Dp-brane throat. For a sufficiently small tempera- 
ture or a sufficiently large separation for the Dq-branes, the probe branes are gravitationally 
attracted towards the horizon but their tension is sufficient to balance this attractive force. 
The probe branes then lie entirely outside of the black hole in what we call a 'Minkowski' 
embedding (see fig. 1). However, above a critical temperature T fun , the gravitational force 

1 Specific examples of this transition were originally seen in [4, 5] and aspects of these transitions in the 
D3/D7 system were independently studied in [6, 7] . Recently, similar holographic transitions have also appeared 
in a slightly different framework [8]. 
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overcomes the tension and the branes are pulled into the horzion. We refer to such configura- 
tions where the branes fall through the horizon as 'black hole' embeddings. In between these 
two phases, there exists a critical solution which just 'touches' the horizon. In [3], we showed 
that in the vicinity of this critical solution the embeddings show a self-similar behaviour. As a 
result, multiple solutions of the embedding equations exist for given temperature in a regime 
close to T fun . Using thermodynamic considerations to select the true ground state then reveals 
a first order phase transition at T fun , where the probe branes jump discontinuously from a 
Minkowski to a black hole embedding. 

In the dual field theory, 2 this phase transition is exemplified by discontinuities in, e.g., 
the quark condensate (ipip) or the contribution of the fundamental matter to the energy 
density. However, the most striking feature of this phase transition is found in the spectrum 
of the mesons, i.e., the quark-antiquark bound states. The latter correspond to excitations 
supported on the probe branes - see, e.g., [9, 10, 11]. In the low-temperature or Minkowski 
phase, the mesons are stable (to leading order within the approximations of large iV c and 
strong coupling) and the spectrum is discrete with a finite mass gap. In the high-temperature 
or black hole phase, stable mesons cease to exist. Rather one finds a continuous and gapless 
spectrum of excitations [12, 13]. Hence the first order phase transition is characterised by the 
dissociation or 'melting' of the mesons. 

This physics is particularly interesting in theories that exhibit a confinement / deconfinement 
phase transition. The dual description of the confining, low-temperature phase involves a 
horizon-free background. At a temperature T deconf the theory undergoes a phase transition at 
which the gluons and the adjoint matter become deconfined, at which point the dual back- 
ground develops a black hole horizon [2]. However, if the mass of the fundamental matter 
is large enough, the branes remain outside the horizon and therefore mesonic bound states 
survive for temperatures T deconf < T < T fun . At T = T fun the branes finally fall into the 
horizon, i.e., the mesons melt. This physics is in qualitative agreement with that observed in 
QCD for heavy-quark mesonic bound states. For example, lattice calculations suggest that 
charmonium states such as the J /ip meson melt at temperatures between 1.65T doconf [14, 15] 
and 2.1T deconf [16], while lattice results for the QCD deconfinement temperature are in the 
range: T dcconf ~ 151 to 192 MeV [17]. Although the holographic description may provide 
some useful geometric intuition for this phenomenon, there are also some caveats that we will 
discuss in due course. 

An overview of the paper is as follows: In section 2, we review the throat geometries for 
black Dp-branes which are dual to (p + l)-dimensional super- Yang-Mills (SYM) at finite tem- 
perature [18]. Section 3 reviews and expands on the self-similar behaviour of the embeddings 
near the critical solution for general Dp/Dq systems, as originally presented in [3]. In the 
subsequent detailed discussion of the thermodynamics, we focus our attention on the D3/D7 
[4] and D4/D6 [5] cases for concreteness. In section 4, we compute the free energy, entropy 
and energy densities, as well as the speed of sound for the D3/D7 system. We also study the 

2 Recall for these supersymmetric field theories, the fundamental matter includes both fermions and scalars, 
which we will refer to collectively as 'quarks'. 
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meson spectrum on the Minkowski embeddings in this section. This spectrum is related to 
the dynamical stability, or lack thereof, of this phase, as we find that tachyonic modes appear 
where thermodynamic considerations indicate that these embeddings are unstable. Section 5 
repeats the salient calculations for the D4/D6 system. Then section 6 concludes with a dis- 
cussion of results. Finally there are several appendices containing various technical details. 
Appendix A provides an analytic description of the D7-brane embeddings at very high and 
very low temperatures. Then appendix B presents some of the details of the calculation of the 
entropy density contributed by the D7-branes. Appendix C discusses the appearance of the 
'swallow tail' form in the plots of the free energy, e.g., fig. 5. Appendix D provides a calcula- 
tion of the constituent quark mass in the low temperature phase of the fundamental matter. 
Finally, appendix E discusses the holographic renormalization of the D4-brane background. 

2. Black Dp-branes 

In this section we briefly review the relevant aspects of the throat geometries and thermo- 
dynamics of black Dp-branes. This will be of use in subsequent sections, in particular, in 
sections 4 and 5, where we specialise to black D3- and D4-brane backgrounds, respectively. 

2.1 Supergravity Background 

The supergravity solution corresponding to the decoupling limit of N c coincident black Dp- 
branes is, in the string frame (see, e.g., [19] and references therein), 

ds 2 = ir 3 (-fdt 2 + dx 2 p ) + m + u 2 dn 2 8 ^ , 

e* = , C i... p = H- 1 , (2.1) 

where H(u) = (L/u) 7 ~ p and /(it) = 1 — {uq/u) 7 ~ p . The horizon lies at u = uq. The length 
scale L is defined in terms of the string coupling constant g B and the string length £ s : 

L 7 - p = g s N c (47^)^ r( V) /4vr . (2.2) 

For the special case p = 3, L is the radius of curvature for the AdSs x S 5 geometry appearing 
in eq. (2.1). 

According to the general gauge/gravity duality of [18], type II string theory in these 
backgrounds is dual to the super- Yang-Mills SU(N C ) gauge theory on the (p + l)-dimensional 
worldvolume of the Dp-branes. For general p 3), the gauge theory is distinguished from 
the conformal case (p = 3) by the fact that the Yang-Mills coupling g YM is dimensionful. The 
holographic dictionary provides 

g YM = 27r 5s (27r4) p - 3 . (2.3) 

Hence there is a power-law running of the dimensionless effective coupling with the energy 
scale U : 

g 2 cH = g 2 YM N c U p - 3 , (2.4) 
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where U = u/a! by virtue of the usual energy/radius correspondence. The absence of confor- 
mal invariance for the general case is manifested in the dual geometry by the radial variation 
of both the string coupling and the spacetime curvature. The supergravity solution (2.1) is 
a trustworthy background provided that both the curvatures and string coupling are small. 
Hence in these general dualities, the supergravity description is limited to an intermediate 
regime of energies in the field theory or of radial distances in the background. This restriction 
is succinctly expressed in terms of the effective coupling (2.4) as [18]: 



Hence the field theory is always strongly coupled where the dual supergravity description is 
valid. 

With the event horizon at u = uq, Hawking radiation appears in the background with a 
temperature fixed by the surface gravity T = k/2tt. This temperature is identified with that 
of the dual {p+ l)-dimensional gauge theory. In the geometry (2.1), the temperature can also 
be determined by demanding regularity of the Euclidean section obtained through the Wick 
rotation t — ► it E . Then t E must be periodically identified with a period (3 where 



In some cases, one periodically identifies some of the Poincare directions x p in order to 
render the theory effectively lower-dimensional at low energies; a prototypical example is that 
of a D4-brane with one compact space direction - see, e.g., [2, 5]. Under these circumstances 
a different background with no black hole may describe the low-temperature physics, and 
a phase transition at T = T dcconf may occur [2]. In the gauge theory this is typically a 
confinement /deconfinement phase transition for the gluonic (or adjoint) degrees of freedom. 
Throughout this paper we assume that T > T doconf , in which case the appropriate gravitational 
background has an event horizon, as in eq. (2.1). 

2.2 Thermodynamics 

Now as alluded to above, with the Wick rotation t — > it E , the Euclidean path integral yields a 
thermal partition function. Further the Euclidean black hole is interpreted as a saddle-point 
in this path integral and so the gravity action evaluated for this classical solution is interpreted 
as the leading contribution to the free energy, i.e., I E = [3F - see, e.g., [20]. Hence to study the 
gauge theory thermodynamics holographically, one needs to evaluate the supergravity action 
I E for the Euclidean version of the above backgrounds (2.1). This suffers from IR (large 
radius) divergences, but these may be regulated by adding appropriate boundary terms to 
the action. These boundary terms were originally found for asymptotically AdS backgrounds, 
such as the black D3-brane, in [21, 22]. As we discuss in appendix E, similar surface terms 
should exist in the general gauge/gravity dualities to complete the holographic description. 
Here we simply comment that for the black D4-brane, which is the relevant background in 



4 



1 < g eB < N c 



(2.5) 




(2.6) 
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section 5, we are guided in the construction of these counterterms by considering the M5- 
brane counterpart in M-theory. In any event, after including the appropriate boundary terms, 
the Euclidean action is finite. 3 Then with F = TI E and standard thermodynamic relations, 
various thermal quantities can be determined. For example, the entropy S and the energy E 
are computed as: 

S = -%i E = F + TS. (2.7) 

For the black D3-brane background, the length scale (2.2) is given by L A = 4irg s N c £*, 
and the free energy is 

F = -^T* = -^N?T\ (2.8) 

where G is the ten-dimensional Newton's constant. In terms of the string length and coupling, 
the latter is given by: 

16vrG = (27r) 7 ^ ffs 2 . (2.9) 
For the black D4-brane geometry we have L 3 = 7rg s N c £g and 

2KV/-9 2 5 tt 2 

F = —ir^ r T & = -^r^T\ (2.io) 

where as usual A = g 2 M N c denotes the 't Hooft coupling. (The reader is referred to appendix 
E for further discussion of this case.) In general, the free energy for a general black Dp-brane 
geometry can be written as [18, 23] 

F ~N 2 a TP +1 g eS {Tf^ , (2.11) 

where 

g 2 s (T) = AT P ~ 3 = g 2 YM N c T p - 3 (2.12) 

is the effective coupling (2.4) evaluated at the temperature scale U = T. In eq. (2.11), iV c 2 
reflects the number of degrees of freedom in the SU(N C ) gauge theory while T p+1 is the 
expected temperature dependence for a (p+ l)-dimensional theory. However, the dependence 
on g eS is a prediction of the holographic framework for the strongly coupled gauge theory. 
Note that for the conformal case (p = 3), but only for this case, this factor is simply unity 
and so the thermodynamic results can compared to those calculated at weak coupling [24]. 

Another quantity that is often studied in the context gauge/gravity duality is the speed 
of sound, e.g., [25, 26, 27, 28, 29]. While this quantity can be inferred from the pole structure 
of certain correlators [25, 26], it can also be derived from the thermal quantities discussed 
above, with 

2 dP _dP fdE\~ 1 S 



dE dT \dT ) c v ' 



3 For the above backgrounds (2.1) describing the gauge theory on flat p- dimensional space, the action still 
contains an IR divergence, namely a factor of the spatial volume V x = J dPx. In the following, we divide all 
extensive thermodynamic quantities by V x so that we are really looking at densities, e.g., eq. (2.11) really gives 
the free energy per unit p-volume. When we refer to contributions from the brane probes, the relevant volume 
factor is instead that of the defect on which the fundamental matter lives, V x = J d d x. 
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Here we have used the fact that for a system without a chemical potential, the pressure and 
free energy density are identical up to a sign, i.e., P = —F. Hence dP/dT = —dF/dT = S. 
Also we use c v to denote the heat capacity (density), i.e., c v = dE/dT. From eqs. (2.11) 
and (2.12), one finds the simple result that for the strongly coupled gauge theory in (p + 1) 
dimensions 

2 = 5-pfl/3 forp=3, 

s 9 -p | 1/5 for p = 4. 1 ; 

We see above that the conformal result Vg = 1/p is only achieved for p = 3 [25, 26], as 
expected. We note, however, that the p = 1 and 4 backgrounds are related through a simple 
chain of dualities to the AdS4 and AdS7 throats of M2- and M5-branes, respectively. Hence 
for these specific cases with v 2 s = 1/2 and 1/5, the speed of sound reflects the conformal 
nature of the holographic theories dual to these M-theory backgrounds [26]. 

3. Criticality, scaling, and phase transitions in Dp/Dq systems 

We now turn to the systems of interest in this paper: Configurations of probe Dq-branes in the 
backgrounds of black Dp-branes. The addition of the probes in the gravitational description 
is dual to the addition of matter in the fundamental representation in the gauge theory [1]. 
This section is mainly a review of [3] that includes some details that were omitted in that 
reference. We describe the embedding of the Dq-brane, study the critical behaviour and 
analyse the nature of the phase transition for general p and q. The latter involves extending 
the Euclidean techniques of the previous section to the worldvolume action of the Dq-brane, 
to study the thermal properties of the fundamental matter. This discussion naturally leads 
to sections 4 and 5, where we provide a detailed analysis of the D3/D7 and D4/D6 brane 
systems. 

3.1 Dp/Dq brane intersections 

Consider a configuration of N c coincident black Dp-branes intersecting N f coincident Dq- 
branes along d spacelike directions. In the limit N f <C iV c the Dq-branes may be treated as a 
probe in the Dp-brane geometry (2.1), wrapping an S n inside the S 8 ~ p . We will assume that 
the Dq-brane also extends along the radial direction, so that q = d+n+1. The corresponding 
gauge theory now contains fundamental matter propagating along a (d+l)-dimensional defect. 
To ensure stability, we will consider Dp/Dq intersections which are supersymmetric at zero 
temperature. Generally this means that we are interested in q = p+4, p+2 or p, as studied in 
[10, 11]. In this case, the Ramond-Ramond field sourced by the Dp-branes does not couple to 
the Dq-brane. For the two cases of special interest here, the D3/D7 and the D4/D6 systems, 
one has n = 3 and n = 2 respectively. If the appropriate direction along the D4-brane is 
compactified, then both cases can effectively be thought of as describing the dynamics of a 
four-dimensional gauge theory with fundamental matter. 
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3.2 Critical behaviour 

To uncover the critical behaviour of the Dp/Dq brane system, we study the behaviour of the 
probe brane near the horizon, following [30] closely - see also [31]. First it is useful to adapt 
the S 8 ~ p metric in (2.1) to the probe brane embedding, and so we write 

dn 2 8 _ p = d6 2 + sin 2 6 dfl 2 n + cos 2 9 dn 2 _ p _ n . (3.1) 

As described above, the Dq-brane wraps the internal S n with radius sin in this line element. 
Now we zoom in on the near horizon geometry with the coordinates 

-=« + ^, Mt) V - < m > 

with T the temperature defined in (2.6). With these coordinates, the event horizon is at 
z = 0. Further y = denotes the axis running orthogonally to the Dq-brane from the Dp- 
branes. Expanding the metric (2.1) to lowest order in z and y gives Rindler space together 
with some spectator directions: 

ds 2 = -{2nT) 2 z 2 dt 2 + dz 2 + dy 2 + y 2 dQ 2 n + dx\ + ■ • • . (3.3) 

The Dq-brane lies at constant values of the omitted coordinates, so these play no role in 
the following. The Dq-brane embedding is specified by a curve (z(a),y(a)) in the (z,y)- 
plane. Since the dilaton approaches a constant near the horizon, up to an overall constant 
the Dq-brane (Euclidean) action is simply the volume of the brane, namely 

^bulk J daVi 2 + V 2 zy n , (3.4) 

where the dot denotes differentiation with respect to a and the reason for the subscript 'bulk' 
will become clear shortly. This is precisely the action considered in ref. [30]. In the gauge 
z = a the equation of motion takes the form 

zyy + (yy - nz)(l + y 2 ) = 0, (3.5) 

while the gauge choice y = a yields 

yz'z + (nzz - y)(l + z 2 ) = 0. (3.6) 

The two types of embeddings described in the introduction for the full background extend 
to this near-horizon geometry (3.3). Hence the solutions again fall into two classes: 'black 
hole' and 'Minkowski' embeddings - see fig. 1. Black hole embeddings are those for which 
the brane falls into the horizon, and may be characterised by yo, the size of the S n there, 
which is also the size of the induced horizon on the Dq-brane worldvolume. The appropriate 
boundary condition is y = 0, y = yo at z = 0. Minkowski embeddings are those for which 
the brane closes off smoothly above the horizon. These are characterised by the distance of 
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closest approach to the horizon, zq, and satisfy the boundary condition z = 0, z = zq at y = 0. 
There is a simple limiting solution for the equations of motion (3.5): y = sjnz. This critical 
solution just touches the horizon at the point y = z = 0, and so it lies between the above two 
classes. Note that this point is a singularity in the induced metric of the Dq-brane. 

The equation of motion (3.5) enjoys a scaling symmetry: If y = f{z) is a solution, then so 
is y = /(liz)/li for any real positive Li. This transformation rescales zq — ► Zo/li for Minkowski 
embeddings, or yo — > yo/fi for black hole embeddings, which implies that all solutions of a 
given type can be generated from any other one by this scaling transformation. 

Consider now a solution very close to the critical one, y(z) = y/nz + 5,{z). Linearising the 
equation of motion (3.5), one finds that for large z the solutions are of the form £(z) = z v± , 
with 

n y/n 2 - 4(n + 1) 
"± = -2± 2 • ^ 

If n < 4, these exponents have non-vanishing imaginary parts, which leads to oscillatory 
behaviour. It appears that one can also get real exponents with n > 5. However, we will 
show below that no such systems are realized in superstring theory. Hence we will only work 
with n < 4 in the following. In this case it is convenient to write the general solution as 

T- 1 r i 

y = y/nz-\ n asin(alogTz) + 6cos(alogTz) , (3.8) 

(Tz) 2 

where a = y / 4(n + 1) — n 2 /2 and a, b are dimensionless constants determined by zq or yo. It 
is easy to show that under the rescaling discussed above, these constants transform as 




cos(alog^) sin(alog/u) 
-sin(alog^) cos(alog^) 




(3.9) 



This result implies that the solutions exhibit discrete self-similarity and yields critical expo- 
nents that characterise the near-critical behaviour. We refer the reader to [30, 31] for details 
but emphasise that this behaviour depends only on the dimension of the sphere. Hence it is 
universal for all Dp/Dq systems (with n < 4). 

Each near-horizon solution gives rise to a global solution when extended over the full 
spacetime (2.1). Each of these embeddings is characterised two constants, which can be read 
off from its asymptotic behaviour and which can be interpreted as the quark mass M q and 
(roughly) the quark condensate (ipip) in the dual field theory - see below. Both of these 
quantities are fixed by zq or yo- As we will see, the values corresponding to the critical 
solution, M* and (^ip)*, give a rough estimate of the point at which a phase transition 
occurs. 

3.2.1 Real scaling exponents? 

From eq. (3.7), we see that the exponents will be real if the dimension of the internal sphere 
wrapped by the Dq-brane is sufficiently large, i.e., if n > 5. This would be interesting because, 
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whereas the oscillatory behaviour for n < 4 leads to a first order phase transition, as we show 
below, real exponents would seem to lead to a second order phase transition. However, we 
will now argue that (under the same assumption to guarantee stability as above) no such 
analysis can be applied for the Dp/Dq systems that actually arise in superstring theory. 

Choosing a value of n, the dimension of the internal sphere, places restrictions on the 
allowed values of both p and q. The internal S n is a subspace of the spherical part of the 
geometry (2.1) and hence we must have p < 8 — n. We have taken a strict inequality here, i.e., 
we do not consider p = 8 — n, because the size of the n-sphere must vary to have nontrivial 
embeddings and so it can not fill the entire internal (8— p)-sphere. Given that p > 0, 4 we 
need only consider n = 5, 6, 7. 

Next, we note that by T-dualising along the p directions common to both sets of branes, 
the brane configuration is reduced to a DO/Dg' intersection, where q' = n + 1 + (p — d). 
Given the previous restriction on n, we must have q' > 6. Now, if we require as above that 
the intersection be supersymmetric at zero temperature (for stability), then we must have 
q' = 8. Hence the only brane configurations of interest are T-dual to the D0/D8 system. 
However, these configurations are those in which string creation arises through the Hanany- 
Witten effect [32]. In particular, as discussed in [33], the background Ramond-Ramond field 
of the Dp-branes will induce a nontrivial worldvolume gauge field on the Dq-brane. While 
this does not rule out the possibility of interesting embeddings and a possible (second order) 
phase transition, it certainly indicates that the present analysis (with no worldvolume gauge 
fields) does not apply to these systems. For this reason, in the remainder of this paper we 
will concentrate on Dp/Dq systems with n < 4. 

3.3 Phase Transitions 

In order to study the global solutions corresponding to the near horizon solutions of the 
previous subsection it is convenient to introduce an isotropic, dimensionless radial coordinate 
p through 

(n 0/ o) ^ = vT^ + ^u 7 -p-u ~ p . (3.10) 

Note that the horizon is at p = 1. Following the discussion in the previous subsection, 5 we 
assume that the Dp/Dq system under consideration is T-dual to the D3/D7 one, in which 
case (p — d) + (n + 1) = 4. Then the Euclidean Dq-brane action density of N { coincident 
Dq-branes in the black Dp-brane background is 

^ = L ip (if 3 (' - 7^) p " (1 ~ v^FW, (3.U) 

4 No black brane geometry exists for a Euclidean D(-l)-brane. 

5 Above, we pointed out that our present analysis does not apply to Dp/Dq systems T-dual to D0/D8- 
branes. Systems T-dual to DO/DO systems would be trivial for the present purposes as n = 0. Hence those 
T-dual to the D0/D4 or D3/D7 system are the only other possibility with a supersymmetric limit. 
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where x = cos9, x = d-x/dp and we have introduced the normalisation constant 

M = JVfTD ^ +ln " . (3.12) 

Here, T Dq = l/(2Tr£ s ) q g s £ s is the Dq-brane tension and Vt n is the volume of a unit n-sphere. 
Up to a numerical constant of O(l), the normalisation factor is found to be 

M ~ N f N c T d g eS (T)^p^ , (3.13) 

where g cB (T) is the effective coupling (2.12) and we have used the standard gauge/gravity 
relations (2.2) and (2.3). 

The equation of motion that follows from (3.11) leads to the large-p behaviour 6 

m c , 
X = - + - + ■■•■ 3.14 
P P 

Holography relates the dimensionless constants m, c to the quark mass and condensate by 7 
M q = ~9 eB (T)^Tm, (3.15) 

(o m) = „ rV (3 . 16) 

47-p 

Here M q is the mass of the fields in the fundamental hypermultiplets, both the fermions tp 
and the scalars q. The operator O m is a supersymmetric version of the quark bilinear, and it 
takes the schematic form 

O m = ^ + q^q + M q qU, (3.17) 

where $ is one of the adjoint scalars. We will loosely refer to its expectation value as the 
'quark condensate'. A detailed discussion of this operator, including a precise definition, can 
be found in appendix A of ref. [35] . 

Eq. (3.15) implies the relation m ( - 5 ~ p ^ 2 = M/T between the dimensionless quantity m, 
the temperature T and the mass scale 



5 — p 
2 

L ) g att (M q ) 



M^f^l ' ~^T. (3.18) 

2 7 -pttL 



Up to numerical factors, this scale is the mass gap in the discrete meson spectrum at tem- 
peratures well below the phase transition [9, 10, 11, 5]. We shall see below that it is also 
the scale of the temperature of the phase transition for the fundamental degrees of freedom, 
T fun ~ M, since the latter takes place at m ~ 1. 

The key observation [31] is that the values (m, c) of a near-critical solution are linearly 
related to the integration constants fixing the corresponding embedding in the near-horizon 



6 Here we assume n > 1. Otherwise the term multiplied by c is logp/p. 
7 Note that the factor of Nt in the second equation was missing in refs. [3, 34]. 
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region. Combining this with the transformation rule (3.9) for the near-horizon constants (a, b) 

— +1 -+i 

and eliminating fi, we deduce that (m — m*)/zQ and (c — c*)/z 2 are periodic functions 
of (a/2ir) log zq with unit period for Minkowski embeddings, and similarly with zq replaced 
by yo for black hole embeddings. This is confirmed by our numerical results, which will be 
discussed in the next sections and are illustrated for the D3/D7 brane system in figure 3. 

The oscillatory behaviour of m and c as functions of zq or yo implies that for a fixed 
value of m near the critical value, several consistent Dq-brane embeddings are possibile with 
different values of c. Alternatively, one finds the quark condensate is not a single-valued 
function of the quark mass. Physically, the preferred solution will be the one that minimises 
the free energy density of the Dq-brane, F = TI Dq . As with the bulk action, the Dq- 
action (3.11) contains large-radius divergences, as can be seen by substituting the asymptotic 
behaviour (3.14) in eq. (3.11). It therefore needs to be regularised and renormalised. We can 
achieve the former by replacing the upper limit of integration by a finite ultraviolet cut-off 
p max . Then in analogy to the holographic renormalisation of the supergravity action [21, 22], 
boundary 'counter-terms' I bound are added to the brane action 7 bu i k , such that the renormalised 
brane energy -/n q — A>uik H~ Abound is then finite as the cut-off is removed, p max > oc [36]. The 
latter method applies directly to asymptotic AdS geometries, but it can be easily extended to 
the D4/D6 system, as discussed below. We expect that a similar procedure can be developed 
for any Dp/Dq system for which there is a consistent gauge/gravity duality. (In any event, 
the brane action can also be regulated by subtracting the free energy of a fiducial embedding.) 
The details for the D3/D7 and D4/D6 cases are discussed in the following sections and the 
results are presented in figures 5 and 12, respectively. In both cases, we see that as the 
temperature is increased, a first order phase transition occurs by discontinuously jumping 
from a Minkowski embedding (point A) to a black hole embedding (point B). We emphasise 
again that this first order transition is a direct consequence of the multi- valued nature of the 
physical quantities brought on by the critical behaviour described in the previous section. It 
may be possible to access this self-similar region by super-cooling the system (although most 
of the other solutions in this region are dynamically unstable - see below). 

It is interesting to ask if the strong coupling results obtained here could in principle be 
compared with a weak coupling calculation. It follows from our analysis that the free energy 
density takes the form F = ft/Tf(m 2 ), where the function / can only depend on even powers 
of m because of the reflection symmetry x —> ~X- The limit m — ► may be equivalently 
regarded as a zero quark mass limit or as a high-temperature limit. In this limit the brane lies 
near the equatorial embedding \ = 0, which slices the horizon in two equal parts. In general 
/(0) is a non-zero numerical constant; in the D3/D7 case, for example, a straightforward 
calculation yields /(0) = —1/2. Given eq. (3.13), we have that at strong coupling the free 
energy density scales as 

F ~ N f N c T d+1 g eS (T) ^ . (3.19) 

The temperature dependence is that expected on dimensional grounds for a <i-dimensional 
defect, and the N { N C dependence follows from large-iV counting rules. However, the depen- 
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dence on the effective 't Hooft coupling indicates that this contribution comes as a strong 
coupling effect, without direct comparison to any weak coupling result. The same is true for 
other thermodynamic quantities such as, for example, the entropy density S = —dF/dT. We 
remind the reader that the background geometry makes the leading contribution to the free 
energy density (2.11), which corresponds to that coming from the gluons and adjoint matter. 
Recall that only for p = 3 is the effective coupling factor absent in eq. (2.11). Only in this 
case the string coupling result differs from that at weak coupling by a mere numerical factor 
of 3/4 [24]. For the fundamental matter, a similar circumstance arises for d = 1, as would 
be realized with the D1/D5, D2/D4 or D3/D3 systems. In these special cases, the strong 
and weak coupling calculations for the fundamental matter could in principle be compared. 
Hence the D3/D3 system is singled out since such a comparison can be made for both the 
adjoint and fundamental sectors. 

4. The D3/D7 system 

Here we will specialise the above discussion to the D3/D7 system. This intersection is sum- 
marised by the array 

0123456789 
D3: x x x x (4.1) 
D7: xxxxxxxx 

Of course, this is an interesting system because both the gluons and the fundamental fields 
in the gauge theory propagate in 3 + 1 dimensions. 

4.1 D7-brane embeddings 

In the D3/D7 brane system with the radial coordinate defined in (3.10), 

(u p) 2 = u 2 + ^x 4 - ut , (4.2) 
the background metric (2.1) becomes 

t 
f 



2 1 fU0P\ 2 

ds = 2 VT ) 



~dt 2 + fdxl 



+ ^[dp 2 + p 2 dn 2 ] , (4.3) 



where 

f(p) = l-^, f(p) = l + ^. (4.4) 

The coordinates {t, x 1 } parametrise the intersection, while {p, ^5} are spherical coordinates 
on the 456789-directions transverse to the D3-branes. As in eq. (3.1), it is useful to adapt 
the metric on the five-sphere to the D7-brane embedding. Since the D7-brane spans the 
4567-directions, we introduce spherical coordinates {r, ^3} in this space and {R, <f\ in the 
89-directions. Denoting by 6 the angle between these two spaces we then have: 

p 2 = r 2 + R 2 , r = psin6>, R = pcos6, (4.5) 
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and 



dp 2 + p 2 dQ\ = dp 2 + p 2 (d6 2 + sin 2 9 dfi§ + cos 2 6 d<f> 2 ) 
= dr 2 + r 2 dn 2 3 + dR 2 + i? 2 # 2 . 



(4.6) 
(4.7) 



Describing the profile in terms of x(p) = cos6>(p) simplifies the analysis - note that 
X = R/p- With this coordinate choice, the induced metric on the D7-brane becomes 



2 1 fuoPV 
dB = 2 (t) 



- f jdt 2 + fdx 2 



L 2 

1 1 7 + 



2-2 



i 2 x 



d / 9 2 + L 2 (l- X 2 )^ 2 , 



(4.8) 



where, as above, x = dx/dp. Since we are studying static embeddings of the probe brane, 
the equation of motion for x(p) can be derived equally well from the Lorentzian or Euclidean 
action. Here we proceed directly to the latter because it is relevant for the thermodynamic 
calculations in the following. The Euclidean D7-brane action density is 



where 



Jdp(l- p 3 (l - X 2 Wl-X 2 + P 2 X 2 , 



= 2vr 2 jV f T D7 ^ = AiV f iV c 3 
4T 32 



(4.9) 



(4.10) 



is the normalisation constant defined in (3.12). Recall from footnote 3 that I bulk denotes a 
density because we have divided out the volume V x . The equation of motion for x{p) ls then 



1 

P 8 



p 5 (i - x 2 )x 



+p° 



1 

P 8 



3x(l ~ X 2 ) + 2p 2 XX 2 

V 1 - x 2 + p 2 x 2 



V 1 - x 2 + p 2 x 2 

which implies that the field x asymptotically approaches zero as 

X 



= 0, 



m c 

- + - + 

P P A 



(4.11) 



(4.12) 



The dimensionless constants m and c are related to the quark mass and condensate through 
eqs. (3.15) and (3.16) with p = 3 and n = 3: 



M = U ° m 

q 



{O m ) = -2 3 / 2 TT 3 e 2 N f T D7 u 3 c=-lVXN t N c T 3 c. 



In this case m = M/T and eq. (3.18) takes the form 

- _ V2(2^ 2 M q ) 2M, 
vrL 2 ^ 



2vr 



(4.13) 
(4.14) 

(4.15) 



where A = g 2 M N c = 2irg s N c is the 't Hooft coupling. In the last equality, we are relating M 
to the meson mass gap in the D3/D7 theory at zero temperature [9]. 
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The equation of motion (4.11) can be recast in terms of the R and r coordinates, related 
to the p and 9 coordinates via (4.5): 



8 r 



d r R 



+ R 2 YJ y/1 + (d r Rf 



r s R 



(r 2 + R?) 



: y/l + (d r R) 2 , (4.16) 



where the embedding of the D7-brane is now specified by R = R(r). Asymptotically, 




Minkowski embedding 
Critical embedding 



Black hole embedding 



2.5 



R(r) = m + — + ... . (4.17) 

In the limits of large and small m we 
were able to find approximate analytic so- 
lutions for the embeddings - see discus- 
sion below and Appendix A. However, for 
arbitrary m we were unable to find an an- 
alytic solution of eq. (4.11) or (4.16) and 
so we resorted to solving these equations 
numerically. It was simplest to solve for 
Minkowski embeddings using the (R, r) co- 
ordinates with equation of motion (4.16) 
while the (x, p) coordinates were best suited Figure 2: Profiles various D7-brane embeddings in a 
to the black hole embeddings. Our ap- D3-brane background in the (R, r)-plane. The circle 
proach was to specify the boundary con- represents the horizon at p = 1. 
ditions at a minimum radius and then nu- 
merically integrate outward. For the black hole embeddings, the following boundary condi- 
tions were specified at the horizon p min = 1: \ = Xo and dx/dp = for < xo < 1- F° r 
Minkowski embeddings, the following boundary conditions were specified at r m j n = (i.e., 
at the axis x = 1) : R = Ro and d r R = for Rq > 1. In order to compute the constants m, c 
corresponding to each choice of boundary conditions at the horizon, we fitted the solutions 
to the asymptotic form (4.12) for x(p) or (4-17) for R(r). A few characteristic profiles are 
shown in fig. 2. 

Recall that, as elucidated in section 3.2, the black hole and Minkowski embeddings are 
separated by a critical solution which just touches the horizon. This critical embedding is 
characterised by certain critical values of the integration constants, m* and c*. For Minkowski 
embeddings near the critical solution, fig. 3 shows plots of (m — m*)/(Ro — l) 5 / 2 and (c — 
c*)/(i?o — l) 5 / 2 versus \/71og(i?o — l)/47r. In this regime, we may relate the boundary value 
to that in the near horizon analysis with Rq — 1 ~ zq. Here our numerical results confirm 
that, near the critical solution, (m — m*) / z^ 2 and (c — c*)/ z^ 2 are both periodic functions of 
\/71og(,zo)/47r with unit period, as discussed above in section 3.2. This oscillatory behaviour 
of m and c as functions of zq (or yo) implies that the quark condensate is not a single- valued 
function of the quark mass and this is clearly visible in our plots of c versus T/M = 1/m, 
displayed in figure 4. By increasing the resolution in these plots, we are able to follow the 
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two families of embeddings spiralling in on the critical solution, the behaviour predicted by 
the near-horizon analysis. Thermodynamic considerations will resolve the observed multi- 
valuedness by determining the physical solution as that which minimizes the free energy 
density of the D7-branes. As discussed in section 3.3, since the physical parameters are 
multi-valued, we can anticipate that there will be a first order phase transition when the 
physical embedding moves from the Minkowski branch to the black hole branch. We will 
proceed to computing the free energy density in the next subsection. The position of the 
resulting phase transition is indicated in the second plot of fig. 4. 

jn - ttl' c - c* 



(R;-!) 5 ' 2 (Ro-l) 5 '" 2 




Figure 3: Quark mass (a) and condensate (b) as a function of the distance to the horizon Rq — 1 for 
D7-brane Minkowski embeddings in a D3-brane background. Note that near the horizon Ro — 1 ~ zq. 

4.2 D7-brane thermodynamics 

Having discussed the embeddings of the D7-brane in the black D3-brane geometry, we proceed 
to compute the free energy, entropy and energy densities associated with the D7-brane, or 
equivalently, the fundamental fields. We start with the Euclidean D7-brane action (4.9). 
Using the asymptotic behaviour (4.12), we see that the action contains a UV divergence, 
since 

J /"Pmax 1 

"Tp-/ dpp^-p m J (4.18) 

** Pmin 

diverges as the regulator is removed, i.e., p max — ► oo. 

This kind of problem is well-known in the context of the AdS / CFT correspondence and 
was first resolved for the gravity action by introducing boundary counter-terms, which depend 
only on the intrinsic geometry of the boundary metric [21, 22]. These ideas can be generalized 
to other fields in an AdS background, such as a scalar [37] - for a review, see [38]. The latter 
formed the basis for the renormalization of probe brane actions in [36], where the brane 
position or profile is treated as a scalar field in an asymptotically AdS geometry. That is, 
one implicitly performs a Kaluza-Klein reduction of the D7 action to five dimensions so that 
it appears to be a complicated nonlinear action for a scalar field \ propagating in a five- 
dimensional (asymptotically) AdS geometry. One then introduces boundary counter-terms 
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Figure 4: Quark condensate c for a D7 in a D3 background versus T/M. The blue dashed (red 
continuous) curves correspond to the Minkowski (black hole) embeddings. The dotted vertical line 
indicates the precise temperature of the phase transition. 

which are local functionals (polynomials) of the scalar field (and boundary geometry) on an 
asymptotic regulator surface. These terms are designed to remove the bulk action divergences 
that arise as the regulator surface is taken off to infinity, as in eq. (4.18). The D3/D7 system is 
explicitly considered in ref. [36], which also introduces a finite counterterm that ensures that 
the brane action vanishes for the supersymmetric embedding of a D7-brane in an extremal 
D3-background, i.e., eq. (2.1) with uq = and p = 3. In the calculation of [36] the D7-brane 
embedding is specified as 9(p), but this is easily converted to a counter-term action for x{p) 
using the obvious coordinate / field redefinition: ^ — 9 = arcsin \ — X + 1/6 X 3 + ' " " • The final 
result is 

= / dtEd3x ^^ " 2x2 + * 4 ) ' (419) 

where this boundary action is evaluated on the asymptotic regulator surface p = p max intro- 
duced above. The boundary metric 7 at p = p max in the (effective) five-dimensional geometry 
is given by 

ds\l) = \ (^) 2 (fc^ dtl + hp^)dxl) (4.20) 
and so ypy = nQ / o max 4 /(/? max )/(p max )/4L 4 . Evaluating the counter-term action (4.19) with an 



asymptotic profile as in eq. (3.14), one finds 



= -\ [(/W 2 - m 2 f ~ 4mc] . (4.21) 

Here we have divided out the volume factor V x - see footnote 3. Comparing eqs. (4.18) 
and (4.21), one sees that the leading divergence proportional to p^ ax cancels in the sum of 
loi = Ajuik + Abound- As a further check, one can consider the supersymmetric limit u — ► 0, in 
which one must work with a rescaled coordinate p = uqP, since the change of variables (4.2) 
is not well defined at uq = 0. In this limit % = u^mjp = \/2 2ir£ 2 M q /p is an exact solution, 
and one can easily verify that for this configuration I D7 = I bulk + I bound = 0. 

In order to produce a finite integral which is more easily evaluated numerically, it is useful 
to incorporate the divergent terms in the boundary action (4.21) into the integral in eq. (4.9) 
using 

Pmax 4 = / dp 4/9 3 + /9 min 4 , 

J Pmin 
/"Pmax 

PmJ = / dp2p + p mh 2 . (4.22) 

J Pmin 

Then the total action may be written as 

^ = G(m) - - [(p min 2 - m 2 ) 2 - Amc] , (4.23) 



M v ; 4 

where G(m) is defined as 



G(m) = dp L 3 (l - 1) (1 - x 2 ) (1 - x 2 + p 2 x 2 ) 1 ' 2 " P 3 + m 2 p 



(4.24) 



Note that the upper bound for the range of integration has been set to infinity, since the 
integral above is finite. 

From these expressions, the free energy density is given by F = Tim. Now using our 
numerical results, the free energy density is shown as a function of the temperature in the 
first two plots in fig. 5. The second of these shows the classic 'swallow tail' form, typically 
associated with a first order phase transition. To our best numerical accuracy, the phase 
transition takes place at T iun /M = 0.7658 (or m = 1.306), where the free energy curves 
for the Minkowski and black hole phases cross. The fact that the transition is first order 
is illustrated by fig. 4, which shows that the quark condensate makes a finite jump at this 
temperature between the points labelled A and B. Similar discontinuities also appear in other 
physical quantities, like the entropy and energy density, as we now calculate. 

Given the free energy density, a standard identity (2.7) yields the entropy density as 



where we have used the expression uq = irL 2 T from eq. (2.6). Evaluating this expression 
requires a straightforward but somewhat lengthy calculation, which we have relegated to 
appendix B. The final result is 

= -4G(m) + (/> min 2 - m 2 ) 2 - 6mc . (4.26) 

Comparing eqs. (4.23) and (4.26), we see that the entropy and free energy densities are simply 
related as 

4F / 2Nmc\ 

s = --f + (4 - 27) 

The first term above can be recognized as the behaviour expected for a conformal system, 
i.e., a system for which F oc T 4 . Hence the second term can be interpreted as summarising 
the deviation from conformal behaviour. We note that, as illustrated in fig. 4, c vanishes in 
both the limits T — ► and T — > oo and so the deviation from conformality is reduced there. 
More precisely, using the results from appendix A we see that c ~ m at high temperature 
and c ~ 1/m 5 at low temperature. Together with (4.10) this implies that the deviation 
from conformality scales as M 2 /T 2 at high temperature. Conformality is also restored at low 
temperatures but only because both S and F/T approach zero more quickly than T 3 . That 
is, S ~ T 7 /M 4 as T -> 0. 

Finally, the thermodynamic identity E = F + TS = T{I m + S) gives the contribution of 
the D7-brane to the energy density: 



L -3G(m) + ? 



, 2 2^2 20 

(Pmin - "1 ) ~ —mC 



(4.28) 



We evaluated both the expressions (4.26) and (4.28) numerically and plotted S and E in 
fig. 5. In both cases, the phase transition is characterised by a finite jump in these quantities, 
as illustrated by the second plot in each case. However, these plots also show that there is 
a large rise in, say, the entropy density in the vicinity of T iun and that the jump associated 
with the phase transition only accounts for roughly 3% of this total increase. 

We close with a few observations about these results. First, recall from (4.10) that 
J\f ~ XN c N t T 3 so that the leading contribution of the D7-branes to all the various thermody- 
namic quantities will be order XN C N { , in comparison to N 2 for the usual bulk gravitational 
contributions. As noted in [3, 34], the factor of A represents a strong coupling enhancement 
over the contribution over a simple free-field estimate for the A^A^ fundamental degrees of 
freedom. We return to this point below in section 6. 

Next, note that in order for the entropy S = —dF/dT to be positive, the free energy F, 
or equivalently the action 7 D7 , must always be a decreasing function of the temperature. This 
means that the apparent 'kinks' in the plot of these quantities versus the temperature are 
true mathematical kinks and not just very rapid turn overs. An analytic proof of this fact is 
given in appendix C. 

Finally, from the plots of the energy density one can immediately read off the qualitative 
behaviour of the specific heat cy = dE/dT. In particular, note that this slope must become 
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negative as the curves spiral around near the critical solution. Hence the corresponding 
embeddings are thermodynamically unstable. Examining the fluctuation spectrum of the 
branes, we will show that a corresponding dynamical instability sets in at precisely the same 
points. One may have thought that these phases near the critical point could be accessed by 
'super-cooling' the system but this instability severely limits the embeddings which can be 
reached with such a process. 




Figure 5: Free-energy, entropy and energy densities for a D7-brane in a D3-brane background; note 
that AT oc T 3 . The blue dashed (red continuous) curves correspond to the Minkowski (black hole) 
embeddings. The dotted vertical line indicates the precise temperature of the phase transition. 
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4.2.1 Thermodynamic expressions for large T/M 

With precisely m = 0, x(p) = is an exact solution. We denote this solution as the equatorial 
embedding, since the D7-brane remains at the maximal S 3 for all values of p. This embedding 
describes the infinite-temperature limit for massive quarks (or massless quarks for any tem- 
perature), i.e., T/M — > oo. For T/M > 1 or m <C 1, approximate analytic solutions for the 
D7-brane profile can be found by perturbing around the equatorial embedding, as discussed 
in appendix A. The final result is given in eq. (A. 4). In the notation of the appendix, the 
integral (4.24) can be expressed as 



G{m) 



00 dx 

T 

-\ + m 2 G 2 
4 



x [ 1 -j ) ( 1 - ^m 2 x 2 + 2x 2 (d x mxf ) - x + mf 



x 



where we have introduced 

r0 ° dx 



Go 



2 



1--) (-h 2 + 2x 2 (d x x) 2 )+l 



0.413893 . (4.29) 



We were only able to evaluate this integral numerically. 

We are now in a position to evaluate the various thermal quantities given by eqs. (4.23), 
(4.26) and (4.28) in this limit. We find 

/ D7 1 ( „ ~ 1\ f M\ 2 1 /M\ 4 



Af . 2+ (-4G 2 -6o-2,^-j +..., 



AfT 2 V 2 / \ T J 

using /9 min = 1 for black hole embeddings and c ~ —0.456947 from eq. (A. 7). In this high 
temperature limit, the quark mass is negligible and so the first term in these expressions could 
be characterised as conformal behaviour. The remaining contributions are small corrections 
indicating a deviation from this simple behaviour generated by the finite quark mass. This 
is essentially the form expected in the high T limit in finite temperature field theory - for 
example, see [39] and the references therein. 

4.2.2 Thermodynamic expressions for small T/M 

Turning to the opposite, low-temperature limit, i.e., T/M <C 1, the D7-branes lie on flat 
embeddings far from the event horizon, i.e., x — Ro/p to leading order. One can calculate 
perturbative improvements to this simple embedding - see appendix A - but it suffices to 
determine the leading thermodynamic behaviour. We find that 



G(m) = [ 
Jo 



oo 

dr 



r :! ( 1- -i ) v i + (O r Bj- + (/•+ RO r B)(nr - ir 



— ^j. (4.30) 
12 m 4 V ; 
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Then using Rq ~ m and c ~ — l/6m 5 , the thermal densities become 

7 /T\ 4 

(4.31) 



1 /T\ 4 S 2 f T\ A E 7/T xl 



M 12 \M) ' 3 VAf / ' NT 12 VM 

Hence these contributions are going rapidly to zero. Note that they still contain the same 
normalization constant (4.10) and so these densities are still proportional to XN f N c . At low 
temperature, one might have expected that the thermodynamics of the fundamental matter 
is dominated by the low lying-mesons, i.e., the lowest energy excitations in the fundamental 
sector, and so that the leading contributions are proportional to iV f 2 , reflecting the number 
of mesonic degrees of freedom. Such contributions to the thermal densities will arise in the 
gravity path integral in evaluating the fluctuation determinant on the D7-brane around the 
classical saddle-point. As indicated by the iV c and A factors, the leading low-temperature 
contributions above come from the interaction of the (deconfined) adjoint fields and the 
fundamental matter. 

4.2.3 Speed of sound 

As mentioned in section 2.2, the speed of sound is another interesting probe of the deconfined 
phase of the strongly coupled gauge theories. In this section, we calculate the effect of 
fundamental matter on the speed of sound. From eq. (2.13), we must evaluate the D7-branes 
contribution to the total entropy density and the specific heat. The first is already given by 
eq. (4.26) and we denote this contribution as S 7 in the following. From eq. (4.28), the energy 
density can be written as E = —3F — 2J\fTmc. Then recalling M oc T 3 from eq. (4.10), the 
D7-brane contribution to the specific heat can be written as 

c v7 = — = 3S 7 -a — (T 4 mc) , (4.32) 

where we have introduced the dimensionless constant a = XN f N c /16. From the black D3 
background, the free energy of the adjoint fields is given in eq. (2.8). It follows then that the 
adjoint contributions to the entropy and specific heat are: 

S 3 = -y^ c 2 T 3 , c v3 = 3S 3 . (4.33) 
Combining all of these results, we can now calculate the speed of sound 



s s 3 + s 7 



c v c v3 + c v7 
S3 + S7 



3S 3 + 3S 7 -ad T (T 4 mc) 

u 1 + ^MT' mc) 



(4.34) 



Note all of our brane calculations are to first-order in an expansion in N f /N c and hence we have 
applied the Taylor expansion in the last line above, reflecting this perturbative framework 
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e.g., c v7 /c v3 -C 1. Now using various expressions above, as well as m 
e = ^jf-, we may write the final result as 



2M, ? /v / AT and 



Svi 



1 

3 



e 
12^ 



1 r^dc\ 
mc H — mT—— 
3 dT 



(4.35) 



This expression indicates that the D7-brane produces a small deviation away from the con- 
formal result, v 2 = g. 

The result of numerically evaluating 5v 2 as a function of the temperature is given in 
fig. 6. We see that 5v 2 is negative. That is, the fundamental matter reduces the speed 
of sound. Following the discussion below eq. (4.27) one finds that 5v 2 ~ T 4 /M 4 at low 
temperature and 5v 2 ~ M 2 /T 2 at high temperature. Thus we see again that the deviation 
from conformal behaviour vanishes for large and small T. We also note that 5v 2 is largest near 
the phase transition, where it makes a discrete jump. Since we are working in a perturbative 
framework, eq. (4.35) is only valid when this deviation is a small perturbation. By assumption 
e oc N f /N c <C 1 and so this is guaranteed provided the last factor in (4.35) is not large. This 
is indeed satisfied for the thermodynamically favoured embeddings, as illustrated in fig. 6. 
Similar deviations have been investigated in [28] for other gauge/gravity dualities. 

In fig. 6, we have also continued 5v 2 on the disfavoured embeddings beyond the phase 
transition and we see that it diverges (towards — oo) at precisely the points where, e.g., the 
energy density curve turns around - see fig. 5. That is, c v7 diverges at these points, so that 
the perturbative derivation of eq. (4.35) breaks down. Hence our perturbative framework 
does not allow us to investigate interesting effects, as seen in [29]. 



^ (3.62 



-0.03 




to. 3 




(a) 



Figure 6: The deviation of the speed of sound from the conformal value, (a) In the limits T — > 
and T — > oo, Sv^ — > 0. (b) The temperature of the phase transition is marked by the dashed vertical 
line. Note there is a finite discontinuity in the speed of sound at the phase transition. If we follow the 
black hole branch (red line) or the Minkowski branch (dotted blue line) past the phase transition, wc 
find that Sv^ diverges. 
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We see from eq. (4.35) that, for massive quarks, the deviation from the conformal result 
is of proportional to N t /N c , as expected from large- AT, counting rules. However, if M q = 
then the result above vanishes, and so 5v 2 s = 0{Nf /N 2 ) at least. Presumably, this additional 
suppression is due to the fact that for massive quarks conformal invariance is broken explicitly 
at the classical level, whereas if M q = it is broken only at the quantum mechanical level by 
the non-vanishing beta function of the theory in the presence of fundamental matter. This is 
proportional to N t /N c , leading to an additional suppression. In the gravitational description 
this is most easily understood at zero temperature. In this case the D3-brane background is 
exactly AdS§ x S 5 , and the isometries of the first factor correspond to the conformal group 
in four dimensions. Adding D7-brane probes with non-zero quark mass breaks the conformal 
isometries, and hence this effect is visible at order N { /N c . Instead, if M q = then the branes' 
world volume is AdS§ x S* 3 , which preserves all the AdS isometries. Hence in this case one 
must go beyond the probe approximation to see the breaking of conformal invariance, i.e., 
beyond 0(N f /N c ). 

4.3 Meson spectrum 

As discussed earlier, introducing the D7-brane probes into the black D3-brane geometry cor- 
responds to adding dynamical quarks into the gauge theory. The resulting theory has a rich 
spectrum of mesons, i.e., quark-antiquark bound states. Since the mesons are dual to open 
strings with both ends on the D7-brane, the mesonic spectrum can be found by computing the 
spectrum of D7-brane fluctuations. For temperatures below the phase transition, T < T fun , 
corresponding to Minkowski embeddings of the D7-branes, we expect the spectrum to exhibit 
a mass gap and be discrete, just as found at T = [9, 10, 11]; this is confirmed by our cal- 
culations below - similar calculations have also appeared recently in [12]. For temperatures 
above the phase transition, corresponding to black hole embeddings, the spectrum will be 
continuous and gapless. Excitations of the fundamental fields in this phase are however char- 
acterised by a discrete spectrum of quasinormal modes, in analogy with [40]. Investigations 
of the black hole phase appear elsewhere [12, 13]. 

4.3.1 Mesons on Minkowski embeddings 

In this section we compute the spectrum of low-lying mesons corresponding to fluctuations 
of the D7-brane in the black D3-brane geometry (4.3). The full meson spectrum would 
include scalar, vector and spinor modes. For simplicity, we will only focus on scalar modes 
corresponding to geometric fluctuations of the brane about the embeddings determined in 
section 4.1. We will work with the (R, r) coordinates introduced in eq. (4.5), in which case 
the background embedding is given by R = R v (r), <p = 0, where the subscript v now indicates 
that this is the 'vacuum' solution. Explicitly, we consider small fluctuations 5R, 5(j) about the 
background embedding: 

R = R v (r) + 5R, <j) = + 5<j). (4.36) 
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The pullback of the bulk metric (4.3) to this embedding is 



2 1 (uop\ 2 
ds = 2 VT ) 



-Ldt? + jdx 2 3 



L 2 

+ — 

P 2 



(1 + R;)dr 2 + r 2 ^ + 2{d a 5R)R v dx a dr 



+- 



L 2 



(d a 5R)(d b 5R)dx a dx b + + 5R) 2 (d a 50)(d b 50)dx a dx b 



where the indices a, 6 run over all D7 worldvolume directions. Using the DBI action, the 
D7-brane Lagrangian density to quadratic order in the fluctuations is 



^ AR v R v d r (5R) 2 A{5R) 2 40R 2 (5R) 2 



-.10 



-,12 



pi°(l+i? 2 ) 

where £q is the Lagrangian density for the vacuum embedding: 



Co = -T^Vh^k 2 (l - 1) . 



(4.37) 



(4.38) 



Here p 2 = r 2 + i? 2 and /i is the determinant of the metric on the S 3 of unit radius. The metric 
g ab in the first line of (4.37) is the induced metric on the D7-brane with the fluctuations set 
to zero: 



ds\ 9 ) = \ 



Cdt 2 + fdx 2 



+ 



L 2 

Pv 



(1 + R 2 v )dr 2 + r 2 dn 2 



(4.39) 



Note that integration by parts and the equation of motion for R v allowed terms linear in 5R 
to be eliminated from the Lagrangian density. The linearised equation of motion is 



da 



L 2 r 3 ffVh 



_P 2 Ji + hl 

for 5R and 
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8Vh 
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IT Ru ■R j v 



da 



r 3 ff^hR^l + ft^sTdaW) 

Pv 



5R 



(4.40) 



for Scj). Summation over the repeated index a is implied. 
We proceed by separation of variables, taking 



V(r) y* 3 (S 3 ) e -*«V k ' x , 5R = K(r) y* 3 (S 3 ) e^e*"* , 



(4.41) 



where y £a (S 3 ) are spherical harmonics on the S 3 , satisfying V| 3 3^ 3 = —^3(^3 + 2)y l ' A . The 
equation of motion for the angular fluctuations becomes 



Or 



r 3 ffR 2 v 

\/i + R 2 v 



d r V 



+ r 3 fR 2 Jl + R? v 
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Co 2 - k 2 - 
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„9 J 



V = 0, (4.42) 
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while for the radial fluctuations we have: 



r 3 ff 

(1 + «2)3/2 



d r K 



+ 



r 3 f 



2 f 2 „ 2 , 2 

7 V [l* u " k 



Ft! 



P 2 , 



4(4 + 2) 



/ 



i + i?? 



ft. 



In these equations, oj and k are dimensionless and are related to their dimensionful counter- 
parts via 



~,2 



a, =a, F = a,7rT =a, 



(4.43) 



and analogously for k. 

We solve these equations using the shooting method. For each choice of the three- 
momentum k, the angular momentum £3, and the embedding R v {r) (corresponding to one 
value of quark mass and chiral condensate) we solve these equations numerically, requiring 
that with r min 0, V(r min ) = r£ in and d r V{r min ) = for the 5<p fluctuations and 

n{r min ) = r% in and d r U{r mm ) = l z r%£ for 5R. Then, as V(r) ~ Ar t:i + Br'^' 2 and 
TZ(r) ~ Cr t} + Dr~ li ~ 2 for some constants A, B, C, D as r — ► 00, we tune J; 2 to find solutions 
which behave as r^ 3-2 asymptotically. 

At finite temperature, the system is no longer Lorentz invariant and so one must consider 
the precise definition for the meson masses. We define the 'rest mass' of the mesons as the 
energy uj with vanishing three-momentum k in the rest-frame of the plasma. 8 Thus, solving 
the equations of motion (4.42) and (4.43) with k = yields the dimensionless constants a) 2 , 
which then give the rest masses through (4.43). 

Plots of the mass spectrum for these modes are given in figs. 7 and 8. Note that in the 
zero-temperature limit, the 5R and 5(f) spectra coincide with those previously calculated for 
the supersymmetric D3 background [9, 10, 11]. In particular, using (4.15), the lightest meson 
in both spectra has a mass squared matching M 2 ap = 47r 2 M 2 ~ 39.5 M 2 . The degeneracy 
between the two different modes arises because supersymmetry is restored at T = and 
both types of fluctuations are part of the same supermultiplet [9, 10, 11]. At finite T, this 
degeneracy between 5R and 5<fi modes is broken. For example, at the phase transition, the 
mass of the lightest meson is roughly 25% and 50% of its zero-temperature value in the 
5R and 5(f) spectra, respectively. The supersymmetric spectrum also showed an unexpected 
degeneracy in that it only depended on the combination n + where n and £3 are the radial 
and angular quantum numbers characterising the individual excitations [9] . Fig. 8 illustrates 
that this degeneracy is broken at finite temperature, where the masses are shown for all the 
modes with n + £3 = 1 and 2. However, this breaking is not large except near the phase 
transition. 



8 Note that this definition differs from [4, 6] which choose M 2 — — k 2 with u — 0. The latter might better 
be interpreted as the low-lying masses of a confining theory in 2+1 dimensions, in analogy to, e.g., [2, 41]. 
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0.7712 0.7714 0.7716 0.771B 



Figure 7: Mass spectrum M 2 = w 2 |fe = o for the 5(f> fluctuations for Minkowski embeddings in the 
D3/D7 system. 

Both figures show that in general the meson masses decrease as the temperature increases. 
As noted above, the thermal shift of the meson rest mass may be of the order of 25 to 50 
percent at the phase transition. This reduction must reflect in part the decrease in the 
constituent quark mass, discussed in appendix D. However, the lowering of the meson masses 
is actually small relative to that seen for the constituent quark mass. As seen in figure 
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0.7712 0.7711 D.7714 0,7715 D.77L6 0.7717 0.7711 

Figure 8: Mass spectrum M 2 — uj 2 \k=o for the SR fluctuations for Minkowski embeddings in the 
D3/D7 system. Note that some of the modes are tachyonic. 

16, at the phase transition, the latter has fallen to only 2% of its T = value. However, 
the thermal shift of the mesons becomes even more dramatic near the critical solution. In 
particular, embeddings with Rq E (1, 1.07) possess tachyonic 5R fluctuations. Note that 
Rq = 1 corresponds to the critical solutions and the phase transition occurs at i?o ~ 1.15, 
i.e., this is the minimum value of Rq for which the thermodynamically preferred embedding 
is of Minkowski type. As discussed above, the embeddings are not unique in the vicinity of 
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the critical solution and so physical quantities spiral in on their critical values. As observed 
at the end of section 4.2, the spiralling of the energy density leads to a negative specific 
heat and indicates an instability. It is satisfying to note in the second plot of fig. 8 that the 
lowest-lying <5i?-mode becomes tachyonic at precisely the point where the first turn-around 
in the spiral occurs (with T). Hence a dynamical instability is appearing in the Minkowski 
embeddings, in precise agreement with the thermodynamic considerations. In fact the second 
lowest-lying 5R-mode becomes tachyonic at the second turn-around and it seems to suggest 
that at the z'th turn of the spiral, the 5 -R- mode with n = i — 1,1% = becomes tachyonic. We 
have found no other evidence of instabilities in other modes. In particular, we have made a 
detailed examination of the spin-one mesons corresponding to excitations of the worldvolume 
gauge field. In this case, the observed behaviour is very similar to that of the pseudoscalar 6(f) 
modes. It is not surprising that a dynamical instability manifests itself in these <5-R-modes, 
since in the region near the critical solution, the nonuniqueness that brings about the phase 
transition arises precisely because the branes have slightly different radial profiles R(r). 

While a dynamical instability set in for the Minkowski embeddings, in agreement with 
the thermodynamic analysis, it is interesting that this point is away from the phase transition. 
In particular, the Minkowski embeddings with Rq £ [1.07,1.15], namely those between the 
point at which the phase transition takes place and the first turn-around, do not exhibit any 
tachyonic modes. Thus these embeddings are presumably meta-stable and might be reached 
through super-cooling. 





(a) 



(b) 



Figure 9: Dispersion relation uo{k) for Minkowski D7-brane embeddings with (a) Rq = 1.20 (to = 
1.32) and (b) Rq = 2.00 (to = 2.00) in a D3-brane background. The solid blue line corresponds to 5<fi 
fluctuations, whereas the red dashed line corresponds to 5R fluctuations. 



We have also made some preliminary investigations of these low-lying mesons moving 
through the thermal plasma and numerical results are shown in figure 9. For non-relativistic 
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motion (small three- momenta) , we expect that the dispersion relation takes the form 

k 2 



u(k) ~ M + 



(4.44) 



where Mq = Mq(T) is the rest mass calculated above and M^ n = M kin (T) is the effective 
kinetic mass for a moving meson. Although M kin (T) is not the same as Mq(T), for low 
temperatures the difference between the two quantities is expected to be small. For example, 
fitting the small-k results for Q for the lowest <5i?-mode at T/M = 0.5 (or Rq = 2) yields 

u) k 2 

— = 6.084 + 0.076^r + • • • . (4.45) 
MM 2 V ' 

Hence in this case, we find M Q /M ~ 6.084 and M kin /M ~ 6.579. Recall that at T = 0, 
we would have M = M kin = M gap = 2irM ~ 6.283M and so both masses have shifted by 
less than 5%. Note that while the rest mass has decreased, the kinetic mass has increased. 
The latter is perhaps counter-intuitive as it indicates it is actually easier to set the meson in 
motion through the plasma than in vacuum. From a gravity perspective, it is perhaps less 
surprising as the Minkowski branes are bending towards the black hole horizon and so these 
fluctuations experience a greater redshift than in the pure AdSs x S 5 background. 

Examining the regime of large three- 
momenta, we find that tu grows linearly 
with k. Naively, one might expect that 
the constant of proportionality should 
be one, i.e., the speed of light. However, 
one finds that 



.6 



0.4 



k 



v m — H — = — h u 
M M 



M 
T 



(4.46) 



0.2 




M 

with v m < 1, as illustrated in fig. 9. 
There our numerical results show that 

for R = 1.2 (m = 1.32), v m ~ 0.353 Figure 10 . The radial profile for SR with Rq = 2 with 
and Mx/M ~ 4.14 for SR and v m ~ various spatial momentum. From top to bottom, the 
0.350 and Mx/M ~ 4.71 for 5</>, while profiles correspond to: k = 0, 4.96, 18.81, 31.4, 39.2, 
for R = 2 (m = 2) v m ~ 0.884 and 62.7, 94.1. 

Mx/M ~ 2.61 for either type of fluctuation. Note that in fig. 9b the dispersion relations oj{k) 
for 5R and 5(j> are nearly coincident for all k because supersymmetry is being restored at 
low temperatures. Our results show that the strongly coupled plasma has a significant effect 
on reducing the maximum velocity of the mesons. This effect is easily understood from the 
perspective of the dual gravity description. The mesonic states have a radial profile which is 
peaked near Rq, the minimum radius of the Minkowski embedding, as illustrated in fig. 10, 
and so we can roughly think of them as excitations propagating along the bottom of the 
D7-brane. At large k, the speed of these signals will be set by the local speed of light 

f(Ro) 



9tt 

9zz 



--Ro 



f(Ro) 



(4.47) 
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The latter gives c ~ 0.349 for Rq = 1.2 and c ~ 0.882 for Rq = 2, both of which closely match 
our results for v m given above. It is interesting that at finite temperature as k increases, the 
radial profiles of the mesonic states seem to become more peaked towards Rq, as illustrated 
in fig. 10. Recall that at T = 0, these profiles are invariant under boosts in the gauge theory 
directions. Finally we note that we did not discover any simple relation between M\ in 
eq. (4.46) and M and M kin in eq. (4.44). 

Note that with the approximations made here, our analysis reveals no dragging forces 
on these low-lying mesons from the thermal bath. We expect that these would only appear 
through string-loop effects, which in particular would include the Hawking radiation of the 
background black hole. This would parallel the similar findings for the drag force experienced 
by large- J mesons composed of heavy quarks [42] and by heavy quarks themselves [43, 44]. 
These large- J mesons also exhibited a maximum velocity similar to the effect discussed above 
[42]. 

5. The D4/D6 system 

We now turn to the D4/D6 system, described by the array 

0123456789 
DA x x x x x (5.1) 
D6 x x x x xxx 

In the decoupling limit, the resulting gauge theory is five-dimensional super- Yang-Mills cou- 
pled to fundamental hypermultiplets confined to a four-dimensional defect. In order to obtain 
a four-dimensional gauge theory at low energies, one may compactify x 4 , the D4-brane di- 
rection orthogonal to the defect, on a circle. If periodic boundary conditions for the adjoint 
fermions are imposed, then supersymmetry is preserved and the four-dimensional theory 
thus obtained is non-confining. In this case the appropriate dual gravitational background 
at any temperature is (2.1) with x periodically identified. Instead, if antiperiodic bound- 
ary conditions for the adjoint fermions are imposed, then supersymmetry is broken and the 
four-dimensional theory exhibits confinement [2] and spontaneous chiral symmetry breaking 
[5]. The holographic description at zero-temperature consists then of D6-brane probes in a 
horizon-free background, whose precise form is not needed here. At a temperature T doconf 
set by the radius of compactification, the theory undergoes a first order phase transition at 
which the gluons and the adjoint matter become deconfined. In the dual description the low- 
temperature background is replaced by (2.1). If T deconf < T fun , the D6-branes remain outside 
the horizon in a Minkowski embedding, and quark-antiquark bound states survive [5] . As T is 
further increased up to T fun a first order phase transition for the fundamental matter occurs. 

Below we study the thermodynamic and dynamical properties of the D6-branes in the 
black D4 background appearing above the deconfinement phase transition. Along the way we 
will have to introduce boundary terms to regulate the D6-brane brane action. 
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5.1 D6-brane embeddings 



As in section 3.3, we begin by transforming to the coordinate system with radial coordinate p 
defined in (3.10), which is better adapted to study the brane embeddings in the background. 
For p = 4, the radial coordinate is then 



and the black D4-brane metric is 
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2 1 fU P\ 3 / 2 

ds = 2 hr) 
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dr + fdxj 
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( L \ 3/2 ulf~^ 



\uopJ 



2V3 



(5.2) 



[dp 2 + P 2 dnl] , (5.3) 



where we now have f(p) = 1 — l/p 3 and f(p) = l + l/p 3 . From eq. (2.6), the temperature is 
given by 

T=^m 112 . (5.4) 



(5.5) 



4vr VL 3 

We also have the holographic relations for the dual five-dimensional gauge theory 



L 3 = irg s NJ 3 s , 



2 



4tt 2 5s 4 , 



where we remind the reader that the Yang-Mills coupling g YM is now dimensionful. 

The D4/D6 intersection is described by the array (5.1). In analogy to the D3/D7 case, 
we introduce spherical coordinates {r, O2} in the 567-directions, and polar coordinates {R, (j)} 
on the 89-plane. Computing boundary terms is also facilitated by introducing an angular 
coordinate between the r and R directions so that we have, as before, 



r 2 + R 2 , r = p sin I 



R = p cos 1 



and 



dp 2 + p 2 dnj = dp 2 + p 2 {d6 2 + sin 2 9 dfll + cos 2 9 dcj) 2 



dr 2 + r 2 d^ll + dR 2 + R 2 d6 2 . 



(5.6) 

(5.7) 
(5.8) 



Following our analysis for the D3/D7 system, we choose coordinates on the brane such that 
asymptotically the metric naturally splits into a product of the form D4-throat x S 2 . We 
describe the embedding of the D6-brane in terms of \{p) = cos 9(p) ~ n °t e then that x = R/p- 
Later, we will have to regulate the Euclidean D6-brane by adding local counter-terms written 
in terms of this 'field.' The induced metric on the D6-brane is then 
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where, as usual, % = dx/dp. The D6-brane action takes the form 



AT 



J dp P 2 (1 - 1) V(i-x 2 )(i-x 2 + p 2 x 2 ), 



dp 2 + p 2 (i- x 2 )dn 2 

(5.9) 
(5.10) 
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where J\f is given by (3.12) with n = 2: 



M = ^N f T D6 u 3 = ^ N f N c g cB (Tf T 3 , 
where g eS (T) 2 = g 2 M N c T. The resulting equation of motion is 



(5.11) 



x 2 x 



V 1 - x 2 + p 2 x 2 

and x asymptotically approaches zero as 
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(5.12) 
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P P 



(5.13) 



with m and c related to the quark mass and condensate via eqs. (3.15) and (3.16) with 
p = 4, n = 2: 
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7.143 



(5.14) 
(5.15) 

(5.16) 



The scale M is again related to the mass gap in the meson spectrum of the D4/D6 system at 
zero temperature. For either background, the latter must be determined numerically. In the 
case of the super symmetric background, one finds [10, 11]: 



™L = 8tt 2 (1.67) 
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W(M q ) 



~ 131.9 



( M q ^ 2 



M 



~ 0.233. 



(5.17) 



One finds essentially the same result for the confining D4 background [5]. The similarity 
of these results is probably a reflection of the underlying supersymmetric structure of the 
five-dimensional gauge theory. In the confining theory, the lowestdying meson is a pseudo- 
Goldstone boson, whose mass is determined by the Gell-Mann-Oakes-Renner relation, and 
the latter linear form extrapolates directly to the supersymmetric result at large M q [5]. 

The equation of motion (5.12) can of course be recast in terms of the R, r coordinates as 



d r 



1 - 



1 



d r R 



pUJ ^l + (d r R) 
which is again suitable to study the Minkowski embeddings. 



6 ^ RyJ 1 + (d r R) 2 , 



(5.18) 
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For arbitrary m we solved for the D6-brane embeddings numerically. Black hole embed- 
dings are most simply described in the x> P coordinates and we used boundary conditions at 
the horizon: \{p = 1) = Xo an d x|p=i = for various < Xo < 1- F° r Minkowski embeddings, 
we used the R, r coordinates and the boundary conditions at the axis were: R(r = 0) = Rq 
and d r R\ r =o = for Rq > 1. We computed m and c by fitting the numerical solutions to 
the asymptotic forms of x an d R given above. In particular, we produced plots of c versus 
T/M, as shown in figure 11. Again by increasing the resolution, we are able to follow the 
two families of embeddings spiralling in on the critical solution. However, thermodynamic 
considerations indicate that a phase transition occurs at the point indicated in the second 
plot. 




T/H T.'K 



Figure 11: Quark condensate c versus temperature T/M for a D6-brane in a D4-brane background. 
The dotted vertical line indicates the precise temperature of the phase transition. 



5.2 D6-brane thermodynamics 

As with the D3/D7 system, we wish to compute the contribution of the fundamental matter 
to the free energy, entropy and energy densities. That is, we will calculate the contributions of 
the D6-brane to the Euclidean path integral. This requires that we regularise and renormalise 
the D6-brane action. We will do this by constructing the appropriate counterterms. 

Using the asymptotic behaviour (5.13) in (5.10) we find that the D6 action contains a 
UV divergence, since 

T C Pmax n 

-jf- ~ / dpp ~ — (5.19) 

diverges for p max — > oo. We expect the counter-terms that must be supplemented to have the 
form J dtvcPxyJdet^ (a + bx 2 + ex 4 )- In the present case, we might expect to pick additional 
factors of e x and e*. In any event, we would choose the constants to eliminate the divergence. 
Further for a supersymmetric embedding, we should be able to construct the counter-term 
action so that the total brane action vanishes. 
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We take as our ansatz for the counter-terms: 



/ bound = 4vrL 3 r D6 [ dt E d 3 x^/fet^e 2a+B *F( X ) , (5.20) 

J P=Pmax 

where B and F{x) are a dimensionless constant and functional of both to be determined. 
We have also defined e 2fJ = ggg; this factor naturally appears in the measure as it is propor- 
tional to the asymptotic volume of the internal S 2 . Now the boundary metric 7^ at p = p max 
in the effective five-dimensional (brane) geometry is given by 

dsHl) = \ (^) 3/2 (g^l dt% + Rp^dx 2 ) (5.21) 
and so ^/det 7 = \ ( n °^ ax ) 3 /(p max )/(p max ) . In this coordinate system we have 



(5.22) 



Now evaluating the counterterm ansatz (5.20) with the supersymmetric background {uq = 0) 
with the profile 9 X = mu o/pi one finds that the leading divergences cancel if B = —2/3 and 
F(0) = —1/3. One also finds that a complete cancellation occurs if we choose 

F{ X ) = - l -(l- X 2 f' 2 . (5.23) 

Thus, the complete counter-term action can be chosen as either of the following: 

/ bound = -^£ 3 T D6 / dt^x^d^ f e 2 °- 2 *l\l - X 2 f 2 , (5.24) 

O J p=Pmax 

/' bound = —~7r L 3 T m [ dt^x^d^e 2 *- 2 ^ (l - 3 - X 2 ) . (5.25) 

o J \ z J p=p max 

In the second expression, we have kept only the terms which contribute to the divergence in 
the small X expansion - the next term of 0(x 4 ) vanishes as p max — ► oo. Computationally, this 
seems like the easier action with which to work; note however that the first form has the nice 
property that, even with finite p max , it produces a precise cancellation for the supersymmetric 
configuration, i.e., I D6 = I hulk + 7 bound = 0. 

Proceeding with /'bound and using (5.13), the boundary term evaluates to 

47rT D6 ( 3 3 2 



/'bound = (^Pmax - -m p max - 3rac J , (5.26) 

where we have divided out the spatial volume V x - see footnote 3. The total action may then 
be written as: 

^ = G{m) - i ^Anm 3 - ^m 2 p min - 3mcj , (5.27) 



9 See the discussion below (4.21). 
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where the integral is defined as 



G{m) = dp \p 2 (l - 1 ) V(1-X 2 )(1-X 2 + P 2 X 2 )-P 2 + \ 

■'Pmin L V P J A . 



(5.28) 



Of course, the free energy follows from this as F = T/ D6 and then one can compute the 
entropy S = —dF/dT and the energy E = F + TS. For the computation of the entropy, 
one must split the free energy into bulk and boundary terms and evaluate the action of the 
derivative on each of the terms, just as was done for the D3/D7 case. We do not present all 
the details of the calculation here but simply give the final result: 



— = -6 G{m) + 2 ( /9 min - -m p min - 4mc 



(5.29) 



where the integral G was defined in (5.28). The contribution of the D6-brane to the energy 
then follows as 
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NT 



21 



5G(m) + - p min - -m p min - —mc 



(5.30) 



Using our numerical results, these thermodynamic quantities are plotted in fig. 12. Again 
the free energy density shows the classic 'swallow tail' form and, to our best numerical accu- 
racy, a first order phase transition takes place at T tun /M = 0.7933 (or m = 1.589), where the 
free energy curves for the Minkowski and black hole phases cross. The fact that the transition 
is first order is illustrated by the entropy and energy densities, which make a finite jump at 
this temperature between the points labelled A and B. 



5.3 Meson spectrum for Minkowski embeddings 

The meson spectrum corresponding to fluctuations of the D6-brane in the black D4-brane 
geometry is computed in the same way as for D3/D7. We focus here on Minkowski embeddings 
for which the spectrum is discrete and stable. The excitations of the black hole embeddings 
will be described by a spectrum of quasinormal modes, as discussed elsewhere [12, 13]. 

We consider small fluctuations 5R, 5(f) about the fiducial embedding, which we now denote 
by Ry, so that the D6-brane embedding is specified by R = R v {r)+5R{x a ) and (f> = 0+5(f)(x a ), 
where R v (r) satisfies (5.18). The pull-back of the bulk metric (5.3) is then 
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(5.31) 



where the metric g is given by 
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Figure 12: Free-energy, entropy and energy densities for a D6-brane in a D4-brane background. 
The blue dashed (red continuous) curves correspond to the Minkowski (black hole) embeddings. The 
dotted vertical line indicates the precise temperature of the phase transition. 



and, as usual, p 2 = r 2 + R 2 . The DBI action yields the D6-brane Lagrangian density to 
quadratic order in the fluctuations 5R, 5(f): 



3R v R v d r {5R) 2 3(5R) 2 2AR 2 {5R) 2 ) 
+ — - : ^ >> 



(d a 5R) 2 
l + Rl 



+ R 2 {d a 



(5.33) 



P 8 V (1 + R 2 V ) Pv 

where h is the determinant of the metric on the S 2 of unit radius, p 2 = r 2 + R 2 , and Cq is 
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the Lagrangian density for the vacuum embedding: 

i 

4 



(5.34) 



Note that terms linear in SR were eliminated from the Lagrangian density C by integration 
by parts and by using the equation of motion (5.18) for R v . Since we are retaining terms only 
to quadratic order in the fluctuations, the metric g v in (5.33) is (5.32) with R = R v and the 
functions / and / in (5.33) and subsequent expressions are evaluated at p v . 
The linearised equations of motion for the fluctuations are then 



da 



ff ( — ^r^^lT^ W 



= 



(5.35) 



for Sip, and 



da 



ff 



L \ 3/2 n2/i/3 r 2^ 



uqPv 



gfd b 5R 



2V3 



l + Rl 



= 6—Vh\/i + m ( i 

Pv 



V 

8-R 2 

Pv 



5R - 6Vhd r 



5R 



for SR. Proceeding via separation of variables, we take 
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where y i2 (S 2 ) are spherical harmonics on the S 2 of unit radius satisfying V| 2 3^ 2 = —^2(^2 + 
l)y £2 . We obtain the radial differential equation 
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and analogously for k. 

We solved (5.38) and (5.37) numerically and determined the eigenvalues Co using a shoot- 
ing method, as was done in the D3/D7 case. The masses are given by M 2 = lo 2 in the 
frame in which the three momentum vanishes: k = 0. The spectra M 2 /M 2 versus T/M 
for the angular fluctuations 5<fi and the radial fluctuations 5R are presented in figs. 13 and 
14, respectively (both for I = 0), and are qualitatively the same as those for the D3/D7 
system: the SR and 5(f) modes become degenerate in the zero-temperature limit, reflecting su- 
persymmetry restoration; in general the meson masses decrease as the temperature increases, 
especially near the critical solution; and the results for 5R fluctuations suggest that a new 
mode becomes tachyonic at each turn-around of the curves. 

6. Discussion 

We have shown that, in a large class of strongly coupled gauge theories with fundamental 
fields, this fundamental matter undergoes a first order phase transition at some high temper- 
ature T fun ~ M, where M is a scale characteristic of the meson physics. As well as giving the 
mass gap in the meson spectrum [9], 1/M is roughly the characteristic size of these bound 
states [45, 11]. In our models, the gluons and other adjoint fields were already in a deconfined 
phase at T fun , so this new transition is not a confinement/deconfinement transition. Neither 
is it a chiral symmetry-restoration phase transition, since the chiral condensate (tpip) oc c that 
breaks the axial U(1) A symmetry does not vanish above T fun . 10 The most striking feature 
of the new phase transition is the change in the meson spectrum and so we refer to it as a 
'dissociation' or 'melting' transition. 

In the low-temperature phase, below the transition, the mesons are deeply bound and 
the spectrum is discrete and gapped. To leading order in the large- -/V c expansion these states 
are absolutely stable, but at higher orders they may decay into other mesons of lower mass or 
glueballs. The leading channel is one-to-two meson decay and after examining the interactions 
in the effective action [9] , we find that parametrically the width of a typical state is given by 
M q /(iVcA 3 / 2 ) ~ M gap /(N C V\). Recall that this is not a confining phase and so we can also 
introduce free quarks into the system. Of course, such a quark is represented by a fundamental 
string stretching between the D7-branes and the horizon. At a figurative level, in this phase, 
we might describe quarks in the adjoint plasma as a 'suspension'. That is, when quarks are 
added to this phase, they retain their individual identities. 

Above the phase transition (i.e., at T > T fun ), the meson spectrum is continuous and gap- 
less. The excitations of the fundamental fields would be characterised by a discrete spectrum 
of quasinormal modes on the black hole embeddings [12, 13]. Investigations of the spectral 
functions [13] show that some interesting structure remains near the phase transition. Some 

10 The large- A^c theories under consideration enjoy an exact U(1)a symmetry, just like QCD at N c = oo. 
However, unlike QCD, they do not possess a non-Abelian SU(Ni)l x SU(N{)n chiral symmetry. Recall also 
that lattice simulations indicate that, in N c = 3 QCD with real-world quark masses, deconfinement and chiral 
symmetry restoration do not occur with a phase transition but through a smooth cross-over [46]. 
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Figure 13: Mass spectrum M 2 — oj 2 \k=o for the 5<fi fluctuations for Minkowski embeddings in the 
D4/D6 system. The dashed vertical line marks the phase transition. 



of these excitations may warrant an interpretation in terms of quasiparticle excitations but in 
any event, there are only a few such states in contrast with the (nominally) infinite spectrum 
of mesons found in the low temperature phase. An appropriate figurative characterisation of 
the quarks in this high temperature phase would be as a 'solution'. If one attempts to inject 
a localised quark charge into the system, it quickly spreads out across the entire plasma and 
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Figure 14: Mass spectrum M 2 = ui 2 \k=o for the 5R fluctuations for Minkowski embeddings in the 
D4/D6 system. The dashed vertical line marks the phase transition while the dotted horizontal line 
marks M 2 = 0. Note that some modes become tachyonic. 

its presence is reduced to diffuse disturbances of the supergravity and worldvolume fields, 
which are soon damped out [12, 13]. 

The physics above is potentially interesting in connection with QCD, since lattice sim- 
ulations indicate that heavy-quark mesons indeed remain bound in a range of temperatures 
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above T deconf . For example, for the lightest charmonium states, the melting temperature may- 
be conservatively estimated to be around 1.65T deconf ~ 249 to 317 MeV [14, 15], depending 
on the precise value of T doconf [17]. Some other studies suggest that the J/ip(lS) state may 
persist to ~ 2.1T deconf ~ 317 to 403 MeV [16]. In the D3/D7 model, we see from fig. 5 
that quark-antiquark bound states melt at T fun ~ 0.766M. The scale M is related to the 
mass M* = M gap of the lightest meson in the theory at zero temperature through eq. (4.15). 
Therefore we have T {un (M*) ~ 0.122M*. For the charmonium states above, taking M* ~ 3000 
MeV gives T fun (cc) ~ 366 MeV. Similarly, for the D4/D6 system we have (5.17) which yields 
M ~ 0.233M*. The transition temperature in this case is then T fun ~ 0.793M ~ 0.186M*, 
which gives T fun (cc) ~ 557 MeV. Hence it is gratifying that these comparisons lead to a 
qualitative agreement with the lattice results. 

Of course, these comparisons must be taken with some caution, since meson bound 
states in Dp/Dq systems are deeply bound, i.e., M* <C 2M q , whereas the binding energy 
of charmonium states is a small fraction of the charm mass, i.e., M c5 ~ 2M C . It might then 
be more appropriate to compare with lattice results for ss bound states which are also seen 
to survive the deconfinement transition. For the 0-meson, whose mass is ~ 1020 MeV, 
the formulas above yield T iun (ss) ~ 124 MeV (D3/D7) and T iun (ss) ~ 188 MeV (D4/D6). 
Lattice simulations suggest that the melting temperature is around 1.4T deconf ~ 211 to 269 
MeV [47, 15]. While again we have qualitative agreement, one must observe that at least for 
the D3/D7 calculation, our result lies below even the lowest estimate for T dcconf ~ 151 MeV. 

An additional caveat is that here we have identified the melting temperature with Tf un , 
above which the discrete meson states disappear. However, the spectral function of some two- 
point meson correlators in the holographic theory still exhibit some broad peaks in a regime 
just above T fun , which suggests that a few bound states persist just above the phase transition 
[13]. This is quite analogous to the lattice approach where similar spectral functions are used 
to examine the existence or otherwise of the bound states. Hence using Tf un above should be 
seen as a (small) underestimate of the melting temperature. 

Before leaving this discussion of comparisons with QCD, we reiterate that the present 
holographic calculations are examining exotic gauge theories and so any agreements above 
must be regarded with a skeptical eye. However, we would also like to point out one simple 
physical parallel between all of these systems. The question of charmonium bound states 
surviving in the quark-gluon plasma was first addressed by comparing the size of the bound 
states to the screening length in the plasma [48]. While the original calculations have seen 
many refinements (see, e.g., [49]), the basic physical reasoning remains sound and so we 
might consider applying the same argument to the holographic gauge theories. Considering 
first the M = 2 SYM theory arising from the D3/D7 system, the size of the mesons can 
be inferred from the structure functions in which the relevant length scale which emerges is 
y/X/M^ [45]. Holographic studies of Wilson lines in a thermal bath [50] show that the relevant 
screening length of the SYM plasma is order 1/T. In fact, the same result emerges from a 
field theoretic scheme of hard-thermal-loop resummation applied to SYM theories [51]. In 
any event, combining these results, the argument that the mesons should dissociate when 
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the screening length is shorter than the size of these bound states yields T ~ M q /\/X Of 
course, the latter matches the results of our detailed calculations in section 4. The same 
reasoning can be applied to the D4/D6 system where the meson size is 0(g eB (M q )/M q ) [11] 
and the screening length is again 0(1 /T) [52]. Hence this line of reasoning again leads to a 
dissociation temperature in agreement with the results of section 5. Therefore we see that 
the same physical reasoning which was used so effectively for the J/ip in the QCD plasma can 
also be used to understand the dissociation of the mesons in the present holographic gauge 
theories. 

One point worth emphasising is that there are two distinct processes that are occurring 
at T ~ M. If we consider, e.g., the entropy density in fig. 5, we see that the phase transition 
occurs in the midst of a cross-over signalled by a rise in 5/T 3 . We may write the contribution 
of the fundamental matter to entropy density as 

S^ = \XN ( N C T* tf(J^) (6.1) 

where H(x) is the function plotted in fig. 5. H rises from 0atx = 0to2asx->oo but 
the most dramatic part of this rise occurs in the vicinity of X = 1. Hence it seems that 
new degrees of freedom, i.e., the fundamental quarks, are becoming 'thermally activated' at 
T ~ M. We might note that the phase transition produces a discontinuous jump in which 
H only increases by about 0.07, i.e., the jump at the phase transition only accounts for 
about 3.5% of the total entropy increase. Thus the phase transition seems to play an small 
role in this cross-over and produces relatively small changes in the thermal properties of the 
fundamental matter, such as the energy and entropy densities. 

As M sets the scale of the mass gap in the meson spectrum, it is tempting to associate 
the cross-over above with the thermal excitation of mesonic degrees of freedom. However, 
the pre-factor \N t N c in (6.1) indicates that this reasoning is incorrect. If mesons provided 
the relevant degrees of freedom, 11 we should have S fun oc iV f 2 . Instead the factor of N t N c 
is naturally interpreted as counting the number of degrees of freedom associated with free 
quarks, with the factor A demonstrating that the contribution of the quarks is enhanced 
at strong coupling. A complementary interpretation of (6.1) comes from reorganizing the 
pre-factor as: 

\N t N c = (gl M N { )N?. (6.2) 

The latter expression makes clear that the result corresponds to the first order correction 
of the adjoint entropy due to loops of fundamental matter. As discussed in [34], we are 
working in a 'not quite' quenched approximation, in that thermal contributions of the D7- 
branes represent the leading order contribution in an expansion in N f /N c , and so fundamental 

11 In fact we will find a contribution proportional to 7V f 2 for the mesons coming from the fluctuation deter- 
minant around the classical D7-brane configuration. One can make an analogy here with the entropy of the 
adjoint fields of N = 4 SYM on S 3 below the deconfinement transition. In this case, the classical gravity saddle- 
point yields zero entropy and one must look at the fluctuation determinant to see the entropy contributed by 
the supergravity modes, i.e., by the gauge-singlet glueballs. 
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loops are suppressed but not completely. In [34], it was shown that the expansion for the 
classical gravitational back-reaction of D7-brane is controlled by XN ( /N C = g 2 M N c . Hence this 
expansion corresponds to precisely the expansion in loops of fundamental matter. However, 
naively the fundamental loops would be suppressed by factors of T 2 /M 2 coming from the 
quark propagators. So from this point of view, the strong coupling enhancement corresponds 
to the fact that such factors only appear as XT 2 /M 2 in eq. (6.1). 

Hence the strongly coupled theory brings together these two otherwise distinct processes. 
That is, at strong coupling, the dissociation of the bound states and the thermal activation 
of the fundamental matter happen at essentially the same temperature. While our discussion 
above focused on the D3/D7 system, the D4/D6 results exhibit the same behaviour. Hence 
this seems to be a universal feature of the holographic gauge theories described by Dp/Dq 
systems. 

The preceding behaviour might be contrasted with that which is expected to occur at 
weak coupling. In this regime, one expects that the melting of the mesons would also be a 
cross-over rather than a (first-order) phase transition. Moreover, the temperature at which 
the mesons dissociate would be T fun ~ E hind ~ gf ff M q . On the other hand, the quarks would 



not be thermally activated until we reach T a< 



2M q , at which point free quark-antiquark 



pairs would be readily produced. Of course, the thermal activation would again correspond 
to a cross-over rather than a phase transition. The key point, which we wish to emphasise, 
is that these two temperatures are widely separated at weak coupling. 

Fig. 15 is an 'artistic' representation of the 
simplest behaviour which would interpolate be- 
tween strong and weak coupling. One might ex- 
pect that the melting point and the thermal ac- 
tivation are very close for g cH 3> 1. The line of 
first order phase transitions must end somewhere 
and so one might expect that it terminates at a 
critical point around g cS ~ 1. Below this point, 
both processes would only represent cross-overs 
and their respective temperatures would diverge 
from one another, approaching the weak coupling 
behaviour described above. 

There are two aspects to enhancement of the 
thermal densities discussed above. First, at strong 
coupling the fundamental matter has a stronger 
effect on the nonabelian plasma than might have 
been otherwise guessed and second the effect is a 
positive one. That is, e.g., the energy and entropy 
densities are raised. Because we are working with 
N { /N c <C 1, the enhancement we observe is a small 
correction to the overall properties of the plasma. In fact, it can be added to a list of such 
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Figure 15: A qualitative representation of 
the simplest possibility interpolating between 
the weak and the strong coupling regimes. 
The solid and the dotted lines correspond 
to T = T fun . At strong coupling this corre- 
sponds to a first-order phase transition (solid 
line), whereas at weak coupling it corre- 
sponds to a cross-over (dotted line). The 
dashed line corresponds to T = T activ . At 
strong-coupling this takes place immediately 
after the phase transition, whereas at weak 
coupling it is widely separated from T fun . 
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correction terms, with others arising as finite-A [53] and finite-A^ effects. Both 12 of these 
types of corrections are expected to raise the entropy and energy densities of the plasma, as 
well. 

As our calculations were also performed in the limit N C ,X — > oo (with N f fixed), it is 
natural to ask how the detailed results of our paper depend on this approximation. First of 
all, the fact that the phase transition is first order implies that it should be stable under small 
perturbations and so its order and other qualitative details should hold within a finite radius 
of the 1/N C ,1/X expansions. Of course, finite- iV c and finite-A corrections may eventually 
modify the behaviour uncovered here. For example, at large but finite iV c the black hole 
will Hawking-radiate and each bit of the brane probe will experience a thermal bath at 
a temperature determined by the local acceleration. This effect becomes more and more 
important as the lower part of a Minkowski brane approaches the horizon, and may potentially 
blur the self-similar, scaling behaviour found here. However, at the phase transition, the 
minimum separation of brane embeddings and the horizon is not parametrically small. For 
example, Rq = 1.1538 at the transition for D3/D7 system. Hence while the Hawking radiation 
can be expected to interfere with the self-similar behaviour near the critical embedding, it 
should not disturb the phase transition for large but finite N c . 

Finite 't Hooft coupling corrections correspond to higher-derivative corrections both to 
the supergravity action and the D-brane action. These may also blur the details of the struc- 
ture discussed above. For example, higher-derivative corrections to the D-brane equation 
of motion are likely to spoil the scaling symmetry of eq. (3.5), and hence the self-similar 
behaviour. These corrections will again become important near the critical solution, for both 
Minkowski and black hole embeddings, since the (intrinsic) curvature of the brane becomes 
large there. However, the phase transition should remain robust for large but finite A because 
at this point, the separation of brane embeddings from the critical solution is not paramet- 
rically small. We illustrated this for the Minkowski embedding at the phase transition of 
D3/D7 above but here we can add the same is true for the black hole embedding at this 
point, which has xo = 0.9427. 

Another significant set of corrections come from the gravitational backreaction of the 
D7-branes (or more generally the probe Dq-branes) on the background spacetime or from 
fundamental loops in the gauge theory. As indicated above, these are dual descriptions of 
the same expansion. Our results only represent the first contribution in an infinite series of 
terms, whose magnitudes are controlled by the ratio N f /N c . Given that low energy QCD has 
N { /N c = 1, it is of particular interest to study holographic theories in Veneziano's limit of 
9ym — ► 0, N c — > oo with both A and N { /N c finite [54]. A variety of attempts have been made 
to construct gravitational backgrounds describing gauge theories in this limit [55]. 

The D2/D6 system provides one interesting background where this limit was studied 
at finite temperature [56]. 13 In particular, it was found that the energy density scales as 

12 At finite N c , the classical black hole would be surrounded by a gas of Hawking radiation which would 
increase both the entropy and energy. 

13 The meson spectrum at T = including the backreaction of the D6-branes has been studied in [57]. 
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F ~ Ny 2 N^ 2 T 3 , which obviously differs from (3.19) with p = 2,d = 2. This discrepancy 
is not at all a contradiction and has the same origin as the discrepancy found for the meson 
spectrum [10, 11]. This is the fact that the calculation in [56] applies in the far infrared of the 
gauge theory, whereas that presented here applies at high temperatures, i.e., at T> g^m- 

We close with a few more observations. Ref. [58] argued for the existence of plasma 
balls in a broad class of confining large- N c theories, which undergo first order deconfinement 
phase transitions. That is, in these theories, one could form metastable, localised lumps 
of deconfined gluon plasma. Their dual description should consist of black holes localised 
along some gauge theory directions. One may imagine an analogous construction for the 
fundamental matter, based on the first order phase transition discussed here. That is, near 
T fun one should be able to construct inhomogenous brane configurations in which only a 
localised region on the branes has fallen through the black hole horizon, i.e., the induced 
brane metric would contain a localised black hole. The dual gauge theory interpretation 
would be in terms of a localised bubble inside of which the fundamental matter has melted. 
Such bubbles may be of interest for understanding how the melting transition actually occurs 
in a dynamical context. 

Finally, we comment on the 'quark condensate' at high temperatures. If one examines 
fig. 4 for example, it is tempting to infer that, since c approaches zero as T — > oo, the 
quark condensate vanishes in this limit. This vanishing would then be in agreement with 
the intuition that at high temperatures the thermal fluctuations should destroy any coherent 
condensate. However, vanishing c is not enough to ensure that (O m ) also vanishes. In fact, if 
we combine eqs. (4.14) and (A. 8), we see that at high temperatures the condensate actually 
grows as 

(O m ) ~ iV c iV f M q T 2 . (6.3) 

At this point, it is important to recall the form of the full operator O m given in eq. (3.17). 
The first two terms are dimension-three operators and so in the high temperature limit we can 
expect the magnitude of typical fluctuations in these to be 0(T 3 ). Further these operators do 
not have a definite sign and so presumably their expectation value vanishes when averaging 
over all fluctuations in the disordered high-temperature system. This, of course, is the basis 
of the intuition that (tpip) — > 0. Now the last term in eq. (3.17) is only a dimension-two 
operator and so we expect thermal fluctuations to be of 0(T 2 ). The key difference in this 
case is that the operator only takes positive real values and so averaging over all fluctuations 
we expect (q'q) oc T 2 . Hence our calculations make a precise prediction for this expectation 
value in the high-temperature phase. Note though that this is a thermal expectation value 
and not a coherent (zero-momentum) condensate, which we expect that we are observing 
with (O m ) ^ at low temperatures. 

Hence it is interesting that the high temperature phase seems to display two distinct 
regimes of behaviour. At very high temperatures, the physics is dominated by incoherent 
thermal fluctuations of the fundamental fields, as expected. However, there is also a regime 
just above the phase transition where the system can support a coherent condensate. This 
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regime would correspond to the region where |c|/T 3 is still growing in fig. 4. Of course, there 
is a cross-over between these two regimes and so there is only a rough boundary. It may 
be natural to define the latter as the point where c is extremized, i.e., T ~ 1.2M. Again, 
this seems to be a universal property of the broad class of holographic theories described by 
Dp/Dq systems. For example, fig. 11 indicates the same behaviour for D6-branes in a D4 
background. 

The above seems to be one more facet of the rich phenomenology which these holographic 
theories display at finite temperature. However, this phenomenology presents several puzzles, 
such as why T fun ~ M rather than M q is the scale at which the bound states melt or at which 
the free quarks are thermally excited. For example, the former seems counterintuitive in view 
of the fact that, in the regime of strong coupling considered here, this temperature is much 
lower than the binding energy of the mesons: 

E bind ~ 2M q - M ~ 2M q . (6.4) 

However, this intuition relies on the expectation that the result of melting a meson is a free 
quark- antiquark pair of mass 2M q . The gravity description makes it clear that this is not the 
case at strong coupling. In fact, the constituent quark mass vanishes when the branes fall 
into the horizon - see appendix D. Rather, in this regime the system is better thought of as 
a strongly coupled liquid of both adjoint and fundamental fields. 

In any event, it is gratifying that the holographic description of these gauge theories 
with fundamental matter provides once more an extremely simple, geometric interpretation of 
some complicated, strong-coupling physics, such as the existence or otherwise of stable quark- 
antiquark bound states above the deconfinement temperature. Other well known examples 
include the geometric characterisations of confinement [2, 59] and chiral symmetry breaking 
[4, 5, 8, 60]. 
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A. Embeddings for high and low temperatures for D3/D7 
A.l High temperatures (black hole embeddings) 

Consider the limit T/M 3> 1. This corresponds to black hole embeddings with m = M/T <C 
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1. As usual, we use the X , p coordinates. Note that the equatorial D7-brane embedding, 
X = 0, is an exact solution of the equation of motion (4.11). To study nearby solutions we 
expand the bulk portion of the D7-brane action (4.9) to quadratic order in X 



i - |y 2 + \e 



thus obtaining the linearised equation of motion: 



On 



1 



P 5 X 



(A.l) 



(A.2) 



To solve this equation, it is useful to make the change of variables x = p 2 so that it becomes: 
x(x 4 - l)(4xx" + 2 X ') + 2x(5x 4 + 3)x + 3(x 4 - l) x = (A.3) 
where X ' = d-x/dx. The solution of this equation satisfying the boundary condition X '\ x =i = 



is 



X 



where F(a, b; c; z) is the hypergeometric function satisfying 

z(l - z)F" + [c - (a + b + l)z]F' - abF = 0. 



(A.4) 



(A.5) 



The overall normalization of the solution is arbitrary since, we are solving a linear equation. 
In the above, we have chosen the normalization such that 



where 



X ~ l/x 1/2 + c/x 3/2 , x -»■ oo 

,_r(^)r(l) 2 



-0.456947. 



(A.6) 



(A.7) 



V2vrr (I) 

The tilde on the solution x and condensate c indicate that these is the solution for unit mass. 
The general solution for arbitrary small mass (or equivalently, high temperatures) is simply 
X = mx and the condensate is given by 



c = mc . 



(A.8) 



A.2 Low temperatures (Minkowski embeddings) 



Low-temperature solutions correspond to Minkowski embeddings in which the D7 probe is 
very far from the horizon: Rq S> 1 or, equivalently, m = M/T> 1. In this case, we expect 
the brane profile to be nearly flat, i.e., R(r) is approximately constant. This motivates the 
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ansatz R(r) = Ro + SR(r), where Rq is a large constant. Substituting into eq. (4.16) and 
expanding to linear order in SR(r) gives: 



r 3 I 1 



1 



(r 2 + R 2 



2\4 



d r (5R) 



r 3 R 



(r 2 + R 2 ) 5 ' 



(A.9) 



Integrating (A.9) and requiring d r (SR)\ r= o = we obtain 



5R(r) 



Ro r , i 

/ ax— 

3 Jo ar 



:! V 1 ( x 2 + #2)4 



i? 2 + 4x 2 



24R§ 



(Ug + * 2 ) 4 i?6 

2(Mj - 1) aretan ( 1 + ^2 + r2) ) + (-1 - 2Hg - 3^) log (l + ^ 



+ (1 - 2i$ + 3i$) log 1 + 



/- \ / 1 + (r 2 + R 2 } 2 ^ 1 



+ 2^ log 



Note that 5-R| r =o = while the limit r — ► oo yields: 

^ (l2^ - 4^) (| - + 12^3 ^ (|^ 



(A.10) 



+^ ( tt^ + tt I log 



8 \3-Rq -Rq 



#o + l 
^-1 



1 1 



+ 



(A.ll) 



Recall that for these embeddings Rq » 1 so indeed 5R <C -Ro- Note that i?(r — ► oo) ~ m+c/r 2 
so that in the large Ro, m limit one has m ~ i?o + 5R|r^oo and c ~ — (6i?o) _1 . For very 
large values of Rq we can expand (A.ll) further to give m ~ i?o + 1/2-Ro as an approximate 
expression for the quark mass. Inverting this relation yields 



Ro c± m 



2m 7 ' 



(A-12) 



We will apply this result in our discussion of the constituent quark mass in appendix D. 



B. Computation of the D7-brane entropy 

In order to evaluate the expression for the D7-brane entropy density, 

dF 2 dF 
S = -df = - nL d^' 



(B.l) 



we split the free energy into a bulk and a boundary contribution. We also write pertinent 
expressions in terms of the dimensionful variables 



p = uq p , c = u c, m = uom , 
to explicitly show the dependences on uq, or, equivalently, the temperature T. 



(B.2) 
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From eq. (4.21), 

Abound = T/ bound = - — T D7 [(p 2 ax " ™ 2 ? ~ 4mc] (B.3) 

so that the boundary contribution to the entropy is 

7r 3 dc 

s hound = - Y L2TmTh d^ > ^ B - 4 ) 

as the quark condensate is the only factor in eq. (B.3) which depends on the position of the 
horizon uq. Note that both of the divergent regulator contributions in eq. (B.3) have been 
eliminated by this differentiation. The bulk contribution to the free energy is given by 

F bulk = T/ bulk = ^T D7 n 4 r ax/M ° dpp* (l X 2 ) (1 - X 2 + P 2 X 2 ) 1 ' 2 ■ (B.5) 

When we differentiate this expression with respect to uq following eq. (B.l), the derivative 
will act in three places: i) the overall factor of uq A ; ii) the explicit (and implicit in p min ) 
appearance of uq in the end-points of the integration; and iii) the field x which is implicitly 
a function of the background mass uq. We consider each of these contributions in turn. First 
one has: 

Si = -2n 3 L 2 T DT u 3 r maX/ "° dpp 3 (l - 1) (1 - x 2 ) (l - X 2 + P 2 X 2 f' 2 ■ (B.6) 

Note that this contribution by itself is divergent in the limit p mM — ► oo. 

Next consider the contributions from the end-points. At the lower end-point, there are 
two possibilities depending on whether the brane ends on the horizon or closes off above the 
horizon. If the brane ends on the horizon, p min = uo and hence this contribution vanishes 
since 3 Uo (/5 min /«o) = 0. (The integrand also vanishes when evaluated at p m in = Anm/^o = 1-) 
If the brane closes off above the horizon, d uo (p m{n /uo) is nonvanishing but this contribution 
vanishes because % = 1 at the end-point. Hence only the upper end-point at p max makes a 
contribution: 



Su = —-tt 3 L 2 T U7 uq 4 



Wl-4)(l-X 2 )(l-X 2 +P 2 X 2 ) 1/2 l x(-M 

V P J Jpmax/UQ V ^0 / 



= ^#T m (p 4 max - m 2 /£J . (B.7) 

where we have substituted the asymptotic expansion (3.14) for \ i n t ne second expression. 

Finally, we consider the contributions from the dependence of x on u o- I n this case, d Ua x 
inside the integral can be considered as a variation <5%. Hence after an integration by parts, 
this derivative yields the bulk equation of motion for \ inside the integral and a boundary 
term coming from the integration by parts. Since x solves the equation of motion, only the 
boundary term contributes to the entropy with 

Pmax/«() 

(B.8) 



~Tr 3 L 2 T D7 u 4 



p / (i - X 2 + P 2 X 2 ) 



- 50 - 



Arguments similar to those above show that the contribution at the lower endpoint vanishes. 
If the brane ends on the horizon, the second factor inside the brackets vanishes and also % 
vanishes at the horizon. If the brane closes off above the horizon, % = 1 at the lower end-point 
and so the numerator in the third factor vanishes and also d uo x = 0. Hence again, only the 
upper end-point contributes to the entropy. In order to correctly evaluate this expression, we 
express the asymptotic expansion (4.12) in terms of m,c: 



X = ^ + ^ + --- (B.9) 
P P A 



Then in eq. (B.8), we have 



fh/uQ Q c/u 3 

X = o p x = o 3 — T~ + " ' 

p z /r 

d UoX = -^-3^ + ^L + ... (B.10) 
p p 6 p 6 

Note that it would be incorrect to evaluate d UQ x — d Uo {rh/p) = because in the integral we 
have assumed that d UQ p = and so d UQ p / 0. Inserting these expansions in (B.8) yields 

S Ui = -^n 3 L 2 T ur [m 2 i5i^ + 6mc-m 4 -u (d Uo c)m] . (B.ll) 

Finally, gathering all the entropy contributions yields: 

S = Si + Sa + Sm + Sbound 



— -- — vr l j d 
2uq 



/ ^ 3 (i-3)(i-x 2 )(i-x 2 + /5 2 x 2 ) 1/2 

J Pmin/U \ P / 

(B.12) 



'Pmin/«0 



Note that the boundary terms have provided precisely the correct p max terms to regulate the 
integral. Hence using (3.12), (4.22) and (4.24), we can express the final result for the entropy 
as 

jf = -4G(m) + (/9 min 2 - m 2 ) 2 - 6mc . (B.13) 
C. Positivity of the entropy 

Here we present an analytic proof that the plot of the Dq-brane probe Euclidean action I Dq 
versus m must exhibit mathematical kinks and not just rapid turn overs. 14 Recall that this is 
necessary for the entropy S = —dF/dT to be positive. We focus here on the case of black hole 
embeddings of the D7-brane in the D3-brane background for concreteness, but the analogous 
arguments applies to Minkowski embeddings and to other Dp/Dq systems. 



14 For the sake of this discussion it is irrelevant whether we plot Iu q versus m or versus 1/m, as in the main 
text. 
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The argument proceeds by thinking of the plot I U7 (m) as a parametric plot (m(xo), Im(xo)), 
where xo> which plays the role of the parameter along the curve, is the value of x at the 'hori- 
zon' p = 1. This is in fact the way we construct the plot: We choose xo as a boundary 
condition at the horizon and we integrate the differential equation 'outwards', thus obtaining 
a solution ip(p;xo), from whose asymptotic behaviour we read off m(xo) and c(xo)- Substi- 
tuting the solution into the D7-brane action we then obtain I D7 (xo)- 

Now the key observation is that if the tangent vector to the curve never vanishes, then 
there can be no kinks. In order to have a kink there must be a point at which both m! = 
dm/dxo and I' = dI D7 /dxo vanish simultaneously. We know that there are certainly an 
infinite number of points at which m' = 0, because close to criticality the function m(xo) is 
an oscillatory function with both maxima and minima. We will now see that at each of these 
points we also have I' D7 = 0. 

The renormalised D7-brane action is I D7 = I bulk + I hound , with 



'bulk 



/ dpjC( X ,x)= dp(l--)p 3 (l-x 2 )Vl-X 2 + P 2 X 2 , (ci) 



and 

4ound = (Pmax 4 - 2m 2 p max 2 - 4?tic + vr£) , (C.2) 

where we have set N = 1 for simplicity. Using the equation of motion, we see that the 
derivative of I D7 is 



j, 5/, 

-*-T\*7 



Ox 







9C dx 
die dxo 



Pmax 



mm 



Pmax 

. (C.3) 



The contribution at p = p min clearly vanishes because /j min = 1 and x(l) = 0. Asymptotically 
we have 

X = - + -,+0{p- A ), (C.4) 
P P A 

and therefore 

*£.-=!: + 4 + 0&,-'). (C.5) 

dxo P P A 
Substituting this into (C.3) we find 

I' D7 = [-mm'p 2 + m 3 m' - 3cm' - mc' + p=f)m ^ . (C.6) 

The derivative of the boundary action is just 

Abound = mm'p^J 2 + mc' + cm - m 3 m' , (C.7) 

so adding everything together we arrive at a simple result in the limit in which the regulator 
is removed: 

/ D7 '^^£I = -2cm'. (C.8) 
dxo 
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This formula is useful for a number of reasons. First, it shows that I' m vanishes if and only 
if m! vanishes, as we wanted to see. Second, applying the chain rule we find 

Physically we expect that c < always, because the brane is attracted to the horizon, and 
this is confirmed by our numerical results. It then follows that J D7 is an increasing function 
of m, or equivalently a decreasing function of 1/m, and hence that the entropy is positive. 
Third, it provides an alternative expression for I D7 , namely 

1 f Xo 

Mxo) = ~2 - 2 7 q dxc(x)m'(x) , (C.IO) 

where we have imposed the boundary condition / D 7(xo = 0) = —1/2, which follows from a 
straightforward calculation of the action of the equatorial embedding. This expression can 
be used to evaluating I m numerically. Moreover, close to criticality one knows the analytic 
form of m(xo) an d c(xo), which should allow one to compute -^(xo) analytically. 

D. Constituent quark mass in the D3/D7 system 

In this section we compute the constituent quark mass M c for temperatures below and near 
the critical temperature for the D3/D7 brane system. A similar analysis has already been 
provided in [43]. 

Our holographic dictionary relates the quark mass M q to the asymptotic constant m with 
eq. (4.13). However this is the bare mass parameter appearing in the microscopic Lagrangian 
of the gauge theory. We must expect the physical or constituent mass of a free quark in the 
deconfined plasma to receive thermal corrections. Since a free quark corresponds to a string 
in the D3-brane geometry hanging from a probe D7-brane (Minkowski embedding) down to 
the horizon, the constituent quark mass corresponds to the energy of this configuration. 

In the notation of the metric (4.3), (4.7), the string worldsheet is extended in the t, R 
directions, localized at r = 0, with induced metric: 

ds 2 = _U^RYf dt 2 + g dR 2 (D1) 



2\LJ f R 2 
The Nambu-Goto string action then becomes 

/ string = -2^2 j dtdRf/fij, (D.2) 

where, since r = 0, / = 1 — 1/i? 4 and / = 1 + 1/i? 4 . Identifying the constituent quark mass 
with minus the action per unit time of this static configuration, we have 



/ 1 \ / 1 N-V2 UQ 




V2 



(D.3) 
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where we recall that i?o = R{? = 0) is the minimal radius reached by the probe brane. Given 
the definition (4.13) for the bare quark mass, we find that 



M c _ 1 



q 



l + i-v/2 



(D.4) 



Plots of M c /M q versus T/M are given in figure 16. In the vicinity of the critical solution, 
there are again multiple embeddings for a fixed value of T/M and so the plots of M c show 
an oscillatory behaviour in this regime. From eq. (D.4), it is clear that as we approach the 
critical solution, i.e., Rq — > 1, the constituent quark mass goes to zero. Note however that 
the phase transition occurs at T/M ~ 0.7658, which corresponds to Rq ~ 1.15 - which is 
marked with the vertical dotted line in fig. 16. Hence the exotic behaviour in the vicinity of 
the critical solution will again not be manifest in the physical system. 

As the temperature goes to zero, M c /M q — ► 1. This is expected, since for small temper- 
atures we have m 3> 1 and we can use the approximate relation (A. 12) in (D.4) to find 

M c y/2 1 5 . . 

1 ~ 1 — — I h ■ ■ ■ fD 5) 

M q m 2m 4 8m 8 v ' ; 

Since m = 1MJ y/\T, this can be finally converted into 

+ •■-. (D.6) 



M c V\T 1 / VAT\ 5 / V\T 

M q ~ ^2M q 2 1 2M q ) 8 I 2M t 



q 



The same expansion appears in [43] but here we have provided an analytic derivation for 
the coefficient of the fourth term, which was obtained in [43] by a numerical fit. Note that 
the two expansions precisely coincide, however, one must replace A — ► 2A above because 
[43] uses a different normalization for the 't Hooft coupling. This difference arises from the 
implicit normalization of the U(N C ) generators: Tr(T a Tfe) = dS a j,. The standard field theory 
convention used in [43] is d = 1/2 while our choice is d = 1, as is prevalent in the D-brane 
literature. For temperatures above the phase transition, the branes fall into the horizon and 
so naively the constituent quark mass vanishes. Rather it is probably inappropriate to speak 
in terms of free quarks in this strongly coupled phase. 



E. Holographic Renormalization of the D4-brane 

Gauge/gravity duality was originally extended to non-AdS backgrounds in [18]. However, 
until recently the discussion of the boundary counter-terms needed for holographic renormal- 
ization [38] was largely limited to asymptotically AdS backgrounds. It was shown that these 
techniques can also be applied in backgrounds describing cascading gauge theory [61]. In prin- 
ciple, we believe it should be possible to extend these techniques to general gauge/gravity 
dualities, in some sense by definition to complete the holographic framework. Here we discuss 
the construction of surface terms which will regulate the Euclidean action of a black Dp-brane 
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Figure 16: Constituent quark mass M c /M q as a function of temperature T/M. The vertical dotted 
line indicates the temperature of the phase transition while the horizontal line indicates that the 
constituent quark mass is roughly M c /M q ~ 0.0212 at the phase transition. Some plots zooming in 
on the spiral behaviour for temperatures slightly above the transition temperature are also shown. 

throat geometry. Again, while formally this may be problematic as generally the supergrav- 
ity description is breaking down in the asymptotic region, some such approach should be 
possible if we believe a gauge/gravity duality exists. We begin with discussion on the D4- 
brane background since this is an interesting place given that it lifts to an (asymptotically) 
AdS7 x S 4 background for which the counter-terms are known. Hence in principle, all we have 
to do is dimensionally reduce the latter to express them in terms of the D4-brane description. 
Given our results for the D4-brane, we make some brief comments on the general Dp-brane 
backgrounds. 

Let us begin by introducing the Euclidean background for a black D4-brane: 



ds 2 



(£) (f(r)dr 2 + dx 2 ) + 



( 



L\ 3 / 2 fdr 2 



r) \f(r) 



_l_ r 




) 



(E.l) 




3 




where 



f(r) = 1 - 



(E.2) 
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and the metric above is in string frame. Recall that the temperature (5.4) and the holographic 
relations (5.5) for the dual five-dimensional gauge theory are given in section 5. 
The string-frame geometry lifts to eleven dimensions as usual: 

ds 2 u = e- 2 *^(ds 2 w ) + e 4 */W , (E.3) 

which for (E.l) yields 

ds 2 = r - (f(r) dr 2 + dx 2 + dz 2 ) + ( -V ^- + L 2 dn 2 
L \r J f(r) 

= (£) 2 (/(«) dr 2 + dx 2 + dz 2 ) + Q + L 2 dtt 2 , (E.4) 

with r/L = (u/L) 2 , f(u) = 1 — (d)/u) 6 , L = 2L and a) 2 = 4Luo- For later discussion, it will 
be convenient to express the throat geometry as 

ds 2 n = e" 2 */ 3 (ds 2 p+2) _ throat + e 2 °L 2 dnl_ p ) + e^dz 2 . (E.5) 

Here the geometry described by ^■5 2 p+ 2)-f/ iroa t replaces the AdS space, while the (p + 2)- 
dimensional field e 2<T describes the running of the internal S 8 ~ p , and L is the scale which we 
will use to replace the AdS scale. 

Now in the standard holographic story for AdS/CFT one refers to the gravity action in 
the same dimension at the AdS space. A similar reduction can be done for the Dp-brane 
throats, i.e., integrate out the internal S 8 ~ p , however it seems more natural to think of them 
as ten-dimensional geometries. Therefore we will consider the bulk action and the Gibbons- 
Hawking surface term are in terms of the full ten (or eleven) dimensions. Note that in the 
usual AdS n x S m examples, these contributions to the action are identical in n and n + m 
dimensions. In particular, note that the sphere factor is constant and so it does not contribute 
to the extrinsic curvature in the Gibbons-Hawking term. So the relevant bosonic terms in 
the Euclidean actions are: 



- - isk / rfI0 ^ (R[S) + 4m)2) - ^i (F<4)>! 



(E.7) 



where we have only kept the terms needed to evaluate the action for the above solution. Also 
we have 16vrGn = (2tt) 8 £ 9 p and 16vrGi = (2vr) 7 ^ 5 2 = (27T) 7 '£%/ 'R n . One subtlety is that 
these two bulk actions are only equal up to an integration by parts. As surface terms play an 
important role in the following, we must keep track of this term. So in reducing the M-theory 
action to the IIA action, one picks up an additional surface term: 

" I d 9 xVh ^e" 2 * n • V<& , (E.8) 

ovrGio J d 
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where h a b denotes the boundary metric in string frame and n is a unit radial vector. Note 
that the norm of the latter is fixed by the ten-dimensional string-frame metric. Now we also 
need the Gibbons-Hawking surface term, which in eleven dimensions is: 



Igh = ~ <fd 10 xVHK n (G) (E.9) 

- j> cfxVhe- 2 * (K 10 (g) - | n ■ . (E.IO) 



2ttR u 



8ttG 

Combining the two ten-dimensional surface terms yields 

I'gh = —tT!- £d 9 xVhe- 2 * (K w {g)+2n-V$>) . (E.ll) 
ovrGio J 

Note that for the D4 throat geometry, the internal S 4 varies with the radial position, and 
so the full ten-dimensional geometry contributes to Kio(g). Hence part of the role of the 
additional term proportional to the radial gradient of <3? is to cancel the sphere contribution, 
as the four-sphere does not contribute in the M-theory calculation. One can check that the 
'unexpected' dilaton term in eq. (E.ll) arises from transforming the standard gravity action 
from Einstein to string frame. 

Now the construction of the remaining boundary counter-terms requires a Kaluza-Klein 
reduction from ten dimensions [62]. For the case of the D4-brane, we can in principle simply 
dimensionally reduce the counter-terms for AdS7, which include a constant or volume term, 
as well as terms proportional to 1Z (the intrinsic curvature) and 1Z 2 . However, we only want 
to consider the D4-brane in Poincare coordinates, i.e., we consider the dual field theory in 
a flat background geometry. Hence the intrinsic curvature contributions will vanish and we 
need only consider the volume term. Note that the prefactor for the AdSy counter-terms 
involves (87rGV) _1 and so we can think that this arose from dimensionally reducing over the 
internal S A . Hence we write the counter-term as: 



let 



— L- f d 4 x^ <f cPxVhI 
87TGn J S i Jd(AdS 7 ) L 

5 OiL 3 



2 8ttG\o 



I d 5 xVhe 4(7 e- 7 */ 3 . (E.12) 
J dM 



So now given the background (E.l), one calculates the Euclidean action Ie as the sum of 
the three terms above in eqs. (E.7), (E.ll) and (E.12). As usual we divide out by the spatial 
volume (see footnote 3), in which case all of the thermodynamic quantities are actually 
densities. In this way we arrive at 

_ n 4 L* /3no 3 _ 2VL»_ 2V 

lE --i^vJ--WGT,¥~~^ r c ' ( } 
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which yields the free energy density given in eq. (2.10). One can also check that this result 
matches that for a planar AdS? black hole [22]. 

Now one can probably extend the counter-term above to general Dp-brane throats. 
The prefactor for the (n — l)-dimensional counter-terms in AdS n x S m examples involves 
(87rG n L) _1 . Hence we have implicitly dimensionally reduced over the internal S m and it 
seems natural that, for the Dp-branes, the prefactor involve Qs-pL 7 ~ p /(8nGio) e^ 8_p ^ cr = 
(87rGp +2 L)- 1 e( 8 - I ') <T . Then it seems the general rule should be that the counter-term takes 
the form 

I ct = A <f d^xVhe^e B \ (E.14) 

where we have written the boundary metric in the string frame, as read off from the ten- 
dimensional or (p + 2)-dimensional string-frame metric, i.e., ds^ p+2 ^_ throat in eq. (E.5). Then 
A and B are dimensionless constants which are chosen experimentally to cancel the relevant 
divergence coming from the bulk and Gibbons-Hawking contributions to the action. 
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